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ABSTRACT 

Corrections are computed to the classical static isotropic solution of general relativity, arising 
from non-perturbative quantum gravity effects. A slow rise of the effective gravitational coupling 
with distance is shown to involve a genuinely non-perturbative scale, closely connected with the 
gravitational vacuum condensate, and thereby, it is argued, related to the observed effective cosmo- 
logical constant. Several analogies between the proposed vacuum condensate picture of quantum 
gravitation, and non-perturbative aspects of vacuum condensation in strongly coupled non-abelian 
gauge theories are developed. In contrast to phenomenological approaches, the underlying functio- 
nal integral formulation of the theory severely constrains possible scenarios for the renormalization 
group evolution of couplings. The expected running of Newton's constant G is compared to known 
vacuum polarization induced effects in QED and QCD. The general analysis is then extended to a 
set of covariant non-local effective field equations, intended to incorporate the full scale dependence 
of G, and examined in the case of the static isotropic metric. The existence of vacuum solutions to 
the effective field equations in general severely restricts the possible values of the scaling exponent 
v. 
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1 Introduction 



Over the last few years evidence has mounted to suggest that quantum gravitation, even though 
plagued by meaningless infinities in standard weak coupling perturbation theory, might actually 
make sense, and lead to a consistent theory at the non-perturbative level. As is often the case 
in physics, the best evidence does not come from often incomplete and partial results in a single 
model, but more appropriately from the level of consistency that various, often quite unrelated, 
field theoretic approaches provide. While it would certainly seem desirable to obtain a closed form 
analytical solution for the euclidean path integral of quantum gravity, experience with other field 
theories suggests that this goal might remain unrealistic in the foreseeable future, and that one 
might have to rely in the interim on partial results and reasoned analogies to obtain a partially 
consistent picture of what the true nature of the ground state of non-perturbative gravity might 
be. 

One aspect of quantum gravitation that has stood out for some time is the rather strident 
contrast between the naive picture one gains from perturbation theory, namely the possibility of 
an infinite set of counterterms, uncontrollable divergences in the vacuum energy of just about any 
field including the graviton itself, and typical curvature scales comparable to the Planck mass [1-3], 
and, on the other hand, the new insights gained from non-perturbative approaches, which avoid 
reliance on an expansion in a small parameter (which does not exist in the case of gravity) and 
which would suggest instead a surprisingly rich phase structure, non-trivial ultraviolet fixed points 
[4-8] and genuinely non-perturbative effects such as the appearance of a gravitational condensate. 
The existence of non-perturbative vacuum condensates does not necessarily invalidate the wide 
range of semi-classical results [9-11] obtained in gravity so far, but re-interprets the gravitational 
background fields as suitable quantum averages, and further adds to the effective gravitational 
Lagrangian the effects of the (finite) scale dependence of the gravitational coupling, in a spirit 
similar to the Euler-Heisenberg corrections to electromagnetism. 

Perhaps the goals that are sometimes set for quantum gravity and related extensions, that is, 
to explain and derive, from first principles, the values of Newton's constant and the cosmological 
constant, are placed unrealistically high. After all, in other well understood quantum field theories 
like QED and QCD the renormalized parameters (a, as, ■ ■■) are fixed by experiments, and no 
really compelling reason exists yet as to why they should take on the actual values observed in 
laboratory experiments. More specifically in the case of gravity, Feynman has given elaborate 
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arguments as to why quantities such as Newton's constant (and therefore the Planck length) might 
have cosmological origin, and therefore unrelated to any known particle physics phenomenon [1]. 

In this paper we will examine a number of issues connected with the renormalization group 
running of gravitational couplings. We will refrain from considering more general frameworks 
(higher derivative couplings, matter fields etc.), and will focus instead on basic aspects of the pure 
gravity theory by itself. Our presentation is heavily influenced by the numerical and analytical 
results from the lattice theory of quantum gravity (LQG), which have, in our opinion, helped 
elucidate numerous details of the non-perturbative phase structure of quantum gravity, and allowed 
a first determination of the scaling dimensions directly in d = 4. The lattice provides a well 
defined ultraviolet regulator, reduces the continuum functional integral to a finite set of convergent 
integrals, and allows statistical field theory methods, including numerical ones, to be used to explore 
the nature of ground state averages and correlations. 

The scope of this paper is therefore to explore the overall consistency of the picture obtained 
from the lattice, by considering a number of core issues, one of which is the analogy with a much 
better understood class of theories, non-abelian gauge theories and QCD (Sec. 2). We will argue 
that, once one takes for granted a set of basic lattice results, it is possible to discuss a number of 
general features without having to explicitly resort to specific aspects of the lattice cutoff or the 
lattice action. For example, it is often sufficient to assume that a cutoff A is operative at very short 
distances, without having to involve in the discussion specific aspects of its implementation. In fact 
the use of continuum language, in spite of its occasional ambiguities when it comes to the proper, 
regulated definition of quantum entities, provides a more transparent language for presenting and 
discussing basic results. 

The second aspect we wish to investigate in this paper is the nature of the rather specific 
predictions about the running of Newton's constant G. A natural starting point is the solution of the 
non-relativistic Poisson equation (Sec. 3), whose solutions for a point source can be investigated for 
various values of the exponent v. We will then show that a scale dependence of G can be consistently 
embedded in a relativistic covariant framework, whose consequences can then be worked out in 
detail for specific choices of metrics (Sec. 4). For the static isotropic metric, we then derive the 
leading quantum correction and show that, unexpectedly, it seems to restrict the possible values for 
the exponent v, in the sense that in some instances no consistent solution to the effective non-local 
field equations can be found unless v~ 1 is an integer. 

To check the overall consistency of the results, a slightly different approach to the solution of the 
static isotropic metric is discussed in Sec. 5, in terms of an effective vacuum density and pressure. 
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Again it appears that unless the exponent v is close to 1/3, a consistent solution cannot be obtained. 
At the end of the paper we add some general comments on two subjects we discussed previously. 
We first make the rather simple observation that a running of Newton's constant will slightly distort 
the gravitational wave spectrum at very long wavelengths (Sec. 6). We then return to the problem 
(Sec. 7) of finding solutions of the effective non-local field equations in a cosmological context [12], 
wherein quantum corrections to the Robertson- Walker metric and the basic Friedman equations 
are worked out, and discuss some of the simplest and more plausible scenarios for the growth (or 
lack thereof) of the coupling at very large distances, past the deSitter horizon. Sec. 8 contains our 
conclusions. 

2 Vacuum Condensate Picture of Quantum Gravitation 

The lattice theory of quantum gravity provides a well defined and regularized framework in 
which non-perturbative quantum aspects can be systematically investigated. 

Let us recall here some of the main results of the lattice quantum gravity (LQG) approach, and 
their relationship to related approaches. 

(i) The theory is formulated via a discretized Feynman functional integral [13-26]. Convergence 
of the euclidean lattice path integral requires in dimensions d > 2 a positive bare cosmological 
constant Ao > [20]. The need for a bare cosmological constant is in line with renormalization 
group results in the continuum, which also imply that radiative corrections will inevitably generate 
a non- vanishing A term. 

(ii) The lattice theory in four dimensions is characterized by two phases, one of which appears for 
G less than some critical value G c , and can be shown to be physically unacceptable as it describes 
a collapsed manifold with dimension d ~ 2. The quantum gravity phase for which G > G c can 
be shown instead to describe smooth four-dimensional manifolds at large distances, and remains 
therefore physically viable. The continuum limit is taken in the standard way, by having the bare 
coupling G approach G c . The two phase structure persists in three dimensions [24], and even at 
d = oo [15], whereas in two dimensions one has only one phase [23]. 

(iii) The presence of two phases in the lattice theory is consistent with the continuum 2 + e 
expansion result, which also predicts the existence of two phases above dimensions d = 2. The 
presence of a nontrivial ultraviolet fixed point in the continuum above d = 2, with nontrivial 
scaling dimensions, relates to the existence of a phase transition in the lattice theory. The lattice 
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results further suggest that the weakly coupled phase is in fact non-perturbatively unstable, with 
the manifold collapsing into a two-dimensional degenerate geometry. The latter phase, if it had 
existed, would have described gravitational screening. 

(iv) One key quantity, the critical exponent v, characterizing the non-analyticity in the vacuum 
condensates at G c , is naturally related to the derivative of the beta function at G c in the 2 + e 
expansion. The value v ~ 1/3 in four dimensions, found by numerical evaluation of the lattice path 
integral, is close but somewhat smaller than the lowest order e expansion result u = \/{d — 2). An 
analysis of the strongly coupled phase of the lattice theory further gives v = at d = oo [15]. 

(v) The genuinely non-perturbative scale £, specific to the strongly coupled phase of gravity 
for which G > G c , can be shown to be related to the vacuum expectation value of the curvature 
via (TV) ~ l/£ 2 , and is therefore presumably macroscopic [34]. It is naturally identified with 
the physical (scaled) cosmological constant A; £ therefore appears to play a role analogous to the 
non-perturbative scaling violation parameter A-^t of QCD. 

(vi) The existence of a non-trivial ultraviolet fixed point (a phase transition in statistical me- 
chanics language) implies a scale dependence for Newton's constant in the physical, strongly coupled 
phase G > G c . To leading order in the vicinity of the fixed point the scale dependence is determined 
by the exponent u, and the overall size of the corrections is set by the condensate scale £. Thus 
in the strongly coupled phase, gravitational vacuum polarization effects should cause the physical 
Newton's constant to grow slowly with distances. 

2.1 Non- Trivial Fixed Point and Scale Dependence of G(/j, 2 ) 

This section will establish basic notation and provide some key results and formulas, some of which 
will be discussed further in the following sections. For more details the reader is referred to the 
more recent papers [13, 15, 12]. 

For the running gravitational coupling one has in the vicinity of the ultraviolet fixed point 



with m = l/£, do > and v c± 1/3 [13]. We have argued that the quantity G c in the above 



expression should in fact be identified with the laboratory scale value, \fG~ c ~ \IG v h ys ~ 1.6 x 

1CT 33 cm, the reason being that the scale £ can be very large. Indeed in [34, 14, 12] it was argued 
that £ should be of the same order as the scaled cosmological constant A. Quantum corrections 
on the r.h.s are therefore quite small as long as k 2 » m 2 , which in real space corresponds to the 
"short distance" regime r <C £. 




(2.1) 
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The above expression diverges as k 2 — > 0, and the infrared divergence needs to be regulated. A 
natural infrared regulator exists in the form of m = l/£, and therefore a properly infrared regulated 
version of the above expression is 

,2 



G(k 2 ) ~ G c 



1 + a 



m 



j_ 

2v 



k 2 + m? 



+ ... 



(2.2) 



with m = l/£ the (tiny) infrared cutoff. Then in the limit of large k 2 (small distances) the correction 
to G{k 2 ) reduces to the expression in Eq. (2.1), namely 

,2\ 2, / 1 m 2 
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+ ... 



whereas its limiting behavior for small k 2 (large distances) is now given by 

k 2 



G(k' 



G c 
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+ ... 



+ 0(r/m 4 ) 



(2.3) 



(2.4) 



,2 i/(l + o ) 

implying that the gravitational coupling approaches the finite value G^ = (1 + a® + . . .) G c , inde- 
pendent of m = l/£, at very large distances r»^. At the other end, for large /c 2 (small distances) 
one has, from either Eqs. (2.2) or (2.1), 

G(k 2 ) ~ G c (2.5) 

k^/m z — » oo 

meaning that the gravitational coupling approaches the ultraviolet (UV) fixed point value G c at 
"short distances" r <C £. Since the theory is formulated with an explicit ultraviolet cutoff, the 
latter must appear somewhere, and indeed G c = A -2 G c , with the UV cutoff of the order of the 
Planck length A^ 1 ~ 1.6 x 10~ 33 cm, and G c a dimensionless number of order one. Note though 
that in Eqs. (2.1) or (2.2) the cutoff does not appear explicitly, it is "absorbed" into the definition 
of G c . 

The non-relativistic, static Newtonian potential is defined as 



(-M)/ 



d 3 k 



e ^x G(k . } __ 



47T 



and therefore proportional to the 3 — d Fourier transform of 
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1 + ao ^ 



+ ... 



(2.6) 



(2.7) 



But, as we mentioned before, proper care has to be exercised in providing a properly infrared 
regulated version of the above expression, which, from Eq. (2.2), reads 

,2 



47T 



47T 



(k 2 + /i 2 ) (k 2 + ^ 2 ) 
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where the limit fi — > should be taken at the end of the calculation. We wish to emphasize here 
that the regulators ji — ► and m are quite distinct. The distinction originates in the condition that 
m arises due to strong infrared effects and renormalization group properties in the quantum regime, 
while n has nothing to do with quantum effects: it is required to make the Fourier transform of the 
classical, Newtonian 4-7r/k 2 well defined. This is an important issue to keep in mind, and to which 
we will return later. 

2.2 Renormalization group properties of G(fj, 2 ) 

This section will discuss the relationship between the running of the coupling G, the renormalization 
group beta function (3(G), the lattice coupling G(A) (the bare coupling, or equivalently, the running 
coupling at the scale of the cutoff A), and the parameter m = l/£. Differentiation of Eq. (2.1) with 
respect to k — ► A gives 

A ^inr = m = - 1 -(g-g c ) + ... (2.9) 

Here and in the following, unless stated otherwise, G will refer to the dimensionless gravitational 
coupling, i.e. G p h ys = A 2 ~ d G(A). In four dimensions, on laboratory scales, yJG v h ys ~ 1CT 33 cm. 
Therefore the exponent \ = — (3'(G C ) is related to the derivative of the beta function for G evaluated 
at the fixed point. Here the dots account for higher order terms not included in either Eq. (2.1) or 
Eq. (2.2). 

The above scaling form for G(k 2 ) and the non-perturbative exponent v are determined as 
follows [34, 13]. Scaling around the fixed point originates in the divergence of the correlation 
length £ = 1/m, related to the bare (lattice) couplings by 

'G(A) - G c 



1 T ! \ ] I T 

m ~ A 

G(A) -» G c 



(2.10) 

where A is the ultraviolet cutoff (the inverse lattice spacing). The continuum limit is approached 
in the standard way by having G — ► G c and A very large, with m kept fixed. This last equation is 
recognized as being just Eq. (2.1) here with the scale k 2 — > A 2 , and solved for the renormalization 
group invariant m (the inverse of the physical correlation length £) in terms of the bare (lattice) 
coupling G(A), at the ultraviolet cutoff scale A. Scaling arguments in the vicinity of the non-trivial 
ultraviolet fixed point then allow one to determine the scaling behavior of correlation functions 
from the critical exponents characterizing the singular behavior of local averages. Since the physical 
quantity m = l/£ is kept fixed and is not supposed to depend on the ultraviolet cutoff A, which is 
sent to infinity, one requires that G(A) change in accordance to 

A^m(A,G(A)) = (2.11) 



It is useful to introduce the dimensionless function F(G) via 

m = r X = AF(G) 



(2.12) 



By differentiation of the renormalization-group invariant quantity m, one then obtains an expression 
for the Callan-Symanzik beta function (3(G) in terms of F. From its definition 

a Ja g{a) = P{G{A)) ' (2 ' 13) 

one has 

« G » - - amho • (2 - 14) 

One concludes that the knowledge of the dependence of m on G, encoded in the function F(G), 
implies a specific form for the (3 function. In terms of the function (3(G) the result of Eq. (2.10) is 
then equivalent to 

13(G) ~ - \ (G - G c ) + 0((G - G c f) (2.15) 

with (3' (G c ) = —\jv. In general one therefore expects the scaling behavior in the vicinity of the 
fixed point 

" ! ~ Aexp (-/ C ^)) G r c „ A l -^- w,G ' , ■ ^ 

where here ~ indicates up to a constant of proportionality. The main conclusion is that the 
function F(G) determines the beta function (3(G), which in turn determines the scale evolution of 
the coupling (obtained from Eq. (2.13), for any //, 

H^-Gin) = (3(G(»)) , (2.17) 

The latter can in principle be integrated in the vicinity of the fixed point, and leads to a definite 
relationship between the relevant coupling G, the renormalization-group invariant (cutoff indepen- 
dent) quantity m = l/£, and the arbitrary sliding scale [i 2 = k 2 , as outlined in the preceding 
section. 

Let us add a comment on the so-called corrections to scaling in the vicinity of the fixed point at 
G c . We note that whereas Eq. (2.2) follows from the result Eq. (2.10) and Eq. (2.9), the infrared 
regulated running coupling of G in Eq. (2.2) is equivalent to assuming the following correction to 
scaling to Eq. (2.10), 



A gfcjgr » (G(K) - GS " 



G(A)^C7 C f\ / G(A) - Gr\ 2l> V a oG c 



/ G(A)-G e \ 
\ a G c J 



Y , 1 (G(k) - G c \ 2v 



(2.18) 
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which gives by differentiation the previously quoted result, Eq. (2.1), plus a small corrections close 
to G c 

(3(G) = -i(G-G c )£ 

where the dots account for higher order terms not included in Eq. (2.2), and, implicitly, in Eq. (2.1) 
as well. 

2.3 Lattice gravity determination of the universal exponent v 

This section summarizes the connection between the lattice regularized quantum gravity path 
integral Z, the singular part of the corresponding free energy F (which determines the scaling 
behavior in the vicinity of the fixed point), and the universal critical exponent v determining the 
scale dependence of the gravitational coupling in Eq. (2.2). For more details the reader is referred 
to [13], and references therein. 

An important alternative to analytic analyses in the continuum is an attempt to solve quantum 
gravity directly via numerical simulations. The underlying idea is to perform the gravitational 
functional integral by discretizing the action on a space-time lattice, and then evaluate the partition 
function Z by summing over a suitable finite set of representative field configurations. In principle 
such a method, given enough configurations and a fine enough lattice, can provide an arbitrarily 
accurate solution to the original quantum gravity theory. 

In practice there are several important factors to consider, which effectively limit the accuracy 
that can be achieved today in a practical calculation. Perhaps the most important one is the 
enormous amounts of computer time that such calculations can use up. This is particularly true 
when correlations of operators at fixed geodesic distance are evaluated. Another practical limitation 
is that one is mostly interested in the behavior of the theory in the vicinity of the critical point 
at G c , where the correlation length £ can be quite large and significant correlations develop both 
between different lattice regions, as well as among representative field configurations, an effect 
known as critical slowing down. Finally, there are processes which are not well suited to a lattice 
study, such as problems with several different length (or energy) scales. In spite of these limitations, 
the progress in lattice field theory has been phenomenal in the last few years, driven in part by 
enormous advances in computer technology, and in part by the development of new techniques 
relevant to the problems of lattice field theories. 

In practice the exponent v in Eqs. (2.10) or (2.18) (and therefore also in Eqs. (2.1) and (2.2), 
which follow from these) is determined from the singularities that arise in the free energy F = 



G - G c 
a G 



2v 



+ 



(2.19) 
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— ylnZ, with the euclidean path integral for pure quantum gravity Z defined as 

Z = J [dg^] e~ A ° I ddx V9+j^Gjd d x y /9R _ 

in the presence of a divergent correlation length £ — > oo in the vicinity of the fixed point at G c . 
This is the scaling hypothesis, the basis of many important results of statistical field theory [38]. 

On purely dimensional grounds, for the singular part of the free energy one has F s i ng (G) ~ 
Standard arguments then give, assuming a divergence £ ~ (G — G c )~ u close to G c , for the first 
derivative of F (here proportional to the average curvature 71(G)) 

"F5G lnZ ~ <Jdx y /g> G~G c AniG - Gc) ■ (2 ' 21) 

An additive constant could appear as well, but the evidence up to now points to this constant being 
zero [13]). Similarly, for the second derivative of F, proportional to the fluctuation in the scalar 
curvature XTl(G), one has 

<UWV>-<IW „ ^(C-CJ-O-* ) . (2.22) 
V dG 2 < f dxyfg > g^g c X1Z v ' y ' 

The above curvature fluctuation is related to the connected scalar curvature correlation at zero 
momentum, 

(Jfe) ^ f dx f dy < ^gRjx) ^gRjy) > c 
< J ' dxyfg> 

and a divergence in the scalar curvature fluctuation is indicative of long range correlations, corre- 
sponding to scale invariance and the presence of massless modes. 

From the computation of such averages one can determine by standard methods the numerical 
values for u, G c and ao to reasonably good accuracy [13]. It is often advantageous to express 
results in the cutoff (lattice) theory in terms of physical (i.e. cutoff independent) quantities. By 
the latter we mean quantities for which the cutoff dependence has been re-absorbed, or restored, 
in the relevant definition. Thus, for example, Eqs. (2.1) and (2.2) will include an overall factor 
of A~ 2 if they refer to the dimensionful, physical Newton's constant; the cutoff is still present, 
but is "hidden" in the definition of physical quantities, and cannot be set equal to infinity as the 
dimensionless fixed point value G c is a finite number. 

As an example, the result equivalent to Eq. (2.21), relating the vacuum expectation value of 
the local scalar curvature (computed therefore for infinitesimal loops) to the physical correlation 
length £ , is 

</y^)> ~ const, fe) 1 -- , (2.24) 
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and which is simply obtained from Eqs. (2.10) and (2.21). Matching of dimensionalities in this last 
equation has been restored by supplying appropriate powers of the Planck length lp = ^G p h ys . 
For v = 1/3 the result of Eq. (2.24) becomes particularly simple [13, 14] 

< fdx x /g R(x) > 1 , 

const. — (2.25) 



< / dxy/g > G->G C lp£ 

The naive estimate, based on a simple dimensional argument, would have suggested the (incorrect) 
result ~ 1/1 p. This shows that v can also play the role of an anomalous dimension, giving the 
magnitude of the deviation from naive dimensional arguments. From the divergences of the free 
energy F one can determine the universal exponent v appearing for example in Eqs. (2.1) and (2.2), 
but not the amplitude ao- The latter requires a direct determination of m = l/£ in terms of bare 
lattice quantities (as in Eq. (2.10)), which we discuss next. 

There are several correlation functions one can compute to extract v and ao directly, either 
through the decay of euclidean invariant correlations at fixed geodesic distance [25], or, equivalently, 
from the correlations of Wilson lines associated with the propagation of heavy spinless particles 
[34]. In either case one expects the scaling result of Eq. (2.10) close to the fixed point, namely 

m ~ A m A\k — k c \ u , A^oo, k — ► k c , m fixed , (2.26) 

with the bare coupling A; (A) = l/(87rG(A)), and A m a calculable numerical constant. Detailed 
knowledge of m{k) allows one to independently estimate the exponent v, but the method is generally 
extremely time consuming (due to the appearance of geodesic distances in the correlation functions), 
and therefore so far not very accurate. 

But, more importantly, from the knowledge of the dimensionless constant A m one can estimate 
from first principles the value of ao in Eqs. (2.1) and (2.2). The first lattice results gave A m ~ 0.56 
[25] and A m ~ 0.87 [34], with some significant uncertainty in both cases (perhaps by as much as 
an order of magnitude, due to the difficulties inherent in computing correlations at fixed geodesic 
distance), which then, combined with the more recent estimate k c ~ 0.0636 and v ~ 0.335 in four 
dimensions [13], gives ao = \/{k c A]t) ~ 42. The rather surprisingly large value for ao appears 
here perhaps as a consequence of the relatively small value of the lattice k c in four dimensions. A 
new determination of ao with significantly reduced errors would clearly be desirable. 

The direct numerical determinations of k c = 1 /8ttG c in d = 3 and d = 4 space-time dimensions 
are in fact quite close to the analytical prediction of the lattice 1/d expansion [15], 

d-2 



kc — 



d^ 



2 d\2 d / 2 
d VdTT 



2/d 

(2.27) 
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The latter gives for a bare A = 1 the estimate k c = \/3/(16 ■ 5 1/4 ) = 0.0724 in d = 4, to be 
compared to the direct determination of k c = 0.0636(11) of [13], and k c = 2 5 / 3 /27 = 0.118 in 
d = 3, to be compared with the direct determination k c = 0.112(5) in [24] (these estimates will be 
compared later in Fig. 1). 

2.4 How many independent bare couplings for pure gravity? 

In this section it will be argued that the scaling behavior of pure quantum gravity is determined by 
one dimensionless combination of Ao and G only. We will then argue that the only sensible scenario, 
from a renormalization group point of view, is one in which the scaled cosmological constant A is 
kept fixed, and only G is allowed to run, as in Eq. (2.2). 

At first it might appear that in pure gravity one has two independent couplings (Ao and G), 
but in reality a simple scaling argument shows that there can only be one, which can be taken 
to be a suitable dimensionless ratio [20]. Consider the (euclidean) Einstein-Hilbert action with a 
cosmological term in d dimensions 

I E [g] = Ao A d J dx Jg- - — L- A d ~ 2 J dx y/g R (2.28) 

Here Ao is the bare cosmological constant and Go the bare Newton's constant, both measured in 
units of the cutoff (we follow customary notation used in cutoff field theories, and denote by A the 
ultraviolet cutoff, not to be confused with the scaled cosmological constant) 3 . Convergence of the 
lattice regulated euclidean path integral requires Ao > [20]. The natural expectation is for the 
bare microscopic, dimensionless couplings to have magnitudes of order one in units of the cutoff, 
Ao ~ Go ~ 0(1). 

Next one rescales the metric so as the obtain a unit coefficient for the cosmological constant 
term, 

4* = r = A - 2/ V' (2.29) 

to obtain 

I E [g] = A d I 'dxy/tf 1 ~^ EI A d - 2 fdx^R' (2.30) 

J 167rG A d J 

The (euclidean) Feynman path integral, defined as 

Z = J [dg^] e-'^l (2.31) 

3 We also deviate in this paper from the convention used in our previous work. Due to ubiquitous ultraviolet 
cutoff A, we reserve here the symbol Ao for the cosmological constant, and A for the scaled cosmological constant 
A ee 8ttG ■ Ao 
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includes as well a functional integration over all metrics, with functional measure given for example 
by [27, 28] 

/ [dg^] = /n [g^)}^ 111 II d 9^(*) (2-32) 

Therefore under a rescaling of the metric the functional measure only picks up a multiplicative 
constant. It does not drop out when computing vacuum expectation values, such as the one in 
Eq. (2.24), but cannot give rise to singularities at finite G such as the ones in Eqs. (2.21) and 
(2.22). 

Equivalently, one can view a rescaling of the metric as simply a redefinition of the ultraviolet 
cutoff A, A -> \l /d A. As a consequence, the non-trivial part of the gravitational functional integral 
over metrics only depends on Ao and Go through the dimensionless combination [20] 

G = G X { o~ 2)/d (2.33) 

The existence of an ultraviolet fixed point is then entirely controlled by this dimensionless parameter 
only, both on the lattice [20, 13, 14] and in the continuum [44] : the non-trivial part of the functional 
integral only depends on this specific combination. One has the Ward identity for the singular part 
of the generating function, 

G »4-^ A "^} F "- (G ' v - ) = (2 ' 34) 

Thus the individual scaling dimensions of the cosmological constant and of the gravitational cou- 
pling constant do not have separate physical meaning; only the relative scaling dimension, as 
expressed through their dimensionless ratio, is physical. 

Physically, the parameter Ao controls the overall scale of the problem (the volume of space-time) , 
while the Go term provides the necessary derivative or coupling term. Since the total volume of 
space-time is normally not considered a physical observable, quantum averages are computed by 
dividing out by the total space-time volume. For example, for the quantum expectation value of 
the Ricci scalar one has the expression of Eq. (2.24). 

Without any loss of generality one can therefore fix the overall scale in terms of the ultraviolet 
cutoff, and set the bare cosmological constant Ao equal to one in units of the ultraviolet cutoff. 4 

The addition of matter field does not change the conclusions of the previous discussion, it is 
just that additional rescalings are needed. Thus for a scalar field with action 

W] = \ J dxj-g (fd^d^ + ml 4? + Rcf) (2.35) 



4 These considerations are not dissimilar from the case of a self-interacting scalar field, where one might want to 
introduce three couplings for the kinetic term, the mass term and the quartic coupling term, respectively. A simple 
rescaling of the field would then reveal that only two coupling ratios are in fact physically relevant. 
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and functional measure (for a single field) 

j[d<f>] = /]J [<?(*)] 1/4 #(*) (2.36) 

the metric rescaling is to be followed by a field rescaling 

<f>'(x) = 4>{x)^' d ' 112 (2-37) 

with the only surviving change being a rescaling of the bare mass tuq — > mo/X^ d . Again the 
scalar functional measure acquires an irrelevant multiplicative factor which does not affect quantum 
averages. 

The same results are obtained if one considers a lattice regularized version of the original 
(euclidean) path integral of Eq. (2.28), which reads [20] 

Z L = J [dl 2 ] e~ lLll2] (2.38) 

with lattice Regge action [16] 

II = A ]T V h (l 2 ) -2k J2 $h(l 2 ) Ml 2 ) (2-39) 
h h 

and regularized lattice functional measure [18, 21, 20] 

[di 2 ] ee / ii4 n [^( s )] ct (^) ( 2 - 4 °) 

ij s 

with ko = l/(167rGo), and 6 a function incorporating the effects of the trian gle inequalities. As 
is customary in lattice field theory, the lattice ultraviolet cutoff is set equal to one (i.e. all lengths 
and masses are measured in units of the cutoff). Convergence of the euclidean lattice functional 
integral requires a positive bare cosmological constant Ao > [20, 22]. On can show again by a 
trivial rescaling of the edges that, as in the continuum, non-trivial part of the lattice regularized 
path integral only depends, in the absence of matter, on the single dimensionless parameter G = 
Gq \ < Q L ~' 2 ^ d . Without loss of generality therefore the bare coupling Ao can be set equal to one 5 . 

The question that remains open is then the following: which coupling should be allowed to 
run within the renormalization group framework? Since the path integral in four dimensions only 



5 In lattice quantum gravity, the average edge length l = (l 2 ) 1 ^ 2 is, for non-singular measures, largely a function 
of Ao. In the large d limit an explicit formula can be given for the measure of Eq. (2.40) with a = 0, 

"2 d\2 d ' 2l2/d 



I 2 - 1 



which agrees well with numerical estimates for finite d [15]. 



(2.41) 
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depends on the dimensionless ratio G 2 = Gq Ao (which is expected to be scale dependent), one has 
several choices; for example G runs and the cosmological constant Ao is fixed. Alternatively, G runs 
and the scaled cosmological constant A = GXq is kept fixed; or G is fixed and A runs etc. 

At first thought, it would seem that the coupling Ao should not be allowed to run, as the overall 
space-time volume should perhaps be considered fixed, not to be rescaled under a renormalization 
group transformation. After all, in the spirit of Wilson [4], a renormalization group transformation 
provides a description of the original physical system in terms of a new coarse-grained Hamiltonian, 
whose new operators are interpreted as describing averages of the original system on the finer scale, 
but within the same physical volume. The new effective Hamiltonian is then supposed to still 
describe the original physical system, but more economically, in terms of a reduced set of degrees 
of freedom. 

These considerations are to some extent implicit in the correct definition of gravitational av- 
erages, for example in Eq. (2.24). Physical, observable averages such as the one in Eq. (2.24) in 
general have some rather non-trivial dependence on the bare coupling Go, more so in the presence 
of an ultraviolet fixed point. Renormalization in the vicinity of the ultraviolet fixed point invariably 
leads to the introduction of a new dynamically generated, non-perturbative scale for G > G c . 

It appears though that the correct answer is that the combination A = 8irG ■ Ao, corresponding 
to the scaled cosmological constant, should be kept fixed, while Newton's constant is allowed to run 
in accordance to the scale dependence obtained from G. The reasons for this choice are three-fold. 
First, in the weak field expansion it is the combination A = GAo that appears as a mass-like term 
(and not Ao or G separately). A similar conclusion is reached if one just compares the appearance of 
the field equations for gravity to say QED (massive via the Higgs mechanism) , or a self-interacting 
scalar field, 

R/iu — \gnvR + Xg^u = 8-kGT^ 
d^F^ + n 2 A u = 4irej u 

d»d^ + m 2 = (2.42) 

Secondly, the scaled cosmological constant represents a measure of physical curvature, as should 
be clear from how the scaled cosmological constant relates for example to the expectation values of 
the scalar curvature at short distances (i.e. for infinitesimally small loops, whose size is comparable 
to the cutoff scale), 

T=rr - R class = 4A (2.43) 

< J dx y/g{x) > 



15 



in the case of pure gravity. But perhaps the most convincing argument for the scaled cosmological 
constant A = 8tt GXo to be kept fixed is given in the following section. 

2.5 The value of £ dilemma - small or large ? 

In this section we will argue that the scale £, which determines the running of G according to 
Eq. (2.2), should be identified with the observed scaled cosmological constant A. 

The lattice quantum gravity result of Eq. (2.24) (and Eq. (2.25) for v = 1/3 ) suggests a 
deep relationship between the correlation length £ = 1/m determining the size of scale dependent 
corrections of Eqs. (2.1) and the curvature. Small averaged curvatures correspond to very large 
length scales £. In gravity, curvature is detected by parallel transporting vectors around closed 
loops. This requires the calculation of a path dependent product of (Lorentz) rotations, R^, 
elements of SO (A) in the euclidean case. On the lattice, the above rotation is directly related to 
the path-ordered (V) exponential of the integral of the lattice affine connection via 



K 



j path T x dx x 

g between simplices 



as 



(2.44) 

p y ' 



as discussed clearly for example in [30, 31], and more recently in [32]. Now, in the strongly coupled 
gravity regime (G > G c ) large fluctuations in the gravitational field at short distances will be 
reflected in large fluctuations of the R matrices, which deep in the strong coupling regime should 
be reasonably well described by a uniform (Haar) measure [15, 29]. 

Borrowing still from the analogy with Yang-Mills theories, one might therefore worry that the 
effects of large strong coupling fluctuations in the R matrices might lead to a phenomenon similar 
to confinement in non-Abelian lattice gauge theories [35] . That this is not the case can be seen 
from the fact that the gravitational analog of the Wilson loop W(F), defined here as a path-ordered 
exponential of the affine connection around a closed planar loop, 



W(T) ~ Tr V exp 



Jc 



(2.45) 



>C 

does not give the static gravitational potential. The static gravitational potential is determined 
instead from the correlation of (exponentials of) geodesic line segments, as in 



cxp 



-M J dryfg^{x) 



dx^ dx v 
dr dr 



(2.46) 



where Mq is the mass of the heavy source, as discussed already in some detail in [33, 34]. Indeed 
a direct lattice calculation of the potential between heavy sources via the correlation of geodesic 
lines showed no sign of confinement [34] . Borrowing from the well-established results in non-abelian 
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lattice gauge theories with compact groups (and to which no exceptions are known), it is easy to 
show that the expected decay of near-planar Wilson loops with area A is then given by 



W(T) ~ exp 



S(C) * 



exp(-A/f) (2.47) 



[36], where A is the minimal physical area spanned by the near-planar loop. The rapid decay of 
the Wilson loop as a function of the area is then seen simply as a general and direct consequence 
of the disorder in the fluctuations of the R matrices at strong coupling. One concludes therefore 
that the Wilson loop in gravity provides a measure of the magnitude of the large-scale, averaged 
curvature, operationally determined by the process of parallel-transporting test vectors around 
very large loops, and which therefore, from the above expression, is computed to be of the order 

r ~ \/i 2 . 

One important assumption in the above result is the identification of the correlation length £ 
in Eq. (2.47) with the correlation length £ in Eqs. (2.10) and (2.24). This is the scaling hypothesis, 
at the basis of most statistical field theory [37, 38, 39] : one assumes that all critical behavior in 
the vicinity of G c is determined by one correlation length, which diverges (in lattice units) at the 
critical point at G c in accordance with Eq. (2.10). As can be seen from Eq. (2.10), the scale £ 
is genuinely non-perturbative, as in non-abelian gauge theories. To determine the actual physical 
value of £ some physical input is needed, as the underlying theory cannot fix it : the ratio of the 
physical Newton's constant to £ 2 can be as small as one desires, provided the bare coupling G is 
very close to its fixed point value G c . 

In conclusion, the above arguments and in particular the result of Eq. (2.47), suggest once more 
the identification of £ with the large scale curvature, the most natural candidate being the (scaled) 
cosmological constant, 

Xphys - ^ OT £ = ~~/T == ( 2 - 48 ) 
S V Aphys 

This relationship, taken at face value, implies a very large, cosmological value for £ ~ 10 28 cm, 
given the present bounds on X p h ys or Hq. Other closely related possibilities may exist, such as 
an identification of £ with the Hubble constant (as measured today) determining the macroscopic 
expansion rate of the universe via the correspondence 

£ - 1/H , (2.49) 

Since this quantity is presumably time-dependent, a possible scenario would be one in which £ _1 = 
Hoo = lim^oo H(t) = \/TT\Ho with = ^P-Xo = ^, where Aq is the observed cosmological 
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constant, and for which the horizon radius is Roo = H^. 

Should Newton's constant run with energy, and if so, according to what law? Newton's constant 
enters the field equations, after multiplication by G itself, as the coefficient of the T^ v term 

R,*, - \g» v R + Xg^ = 8nGT^ (2.50) 

In line with the previous discussion, the running of the (scaled) cosmological term is according to 
the rule A — > l/£, 2 , i.e. no scale dependence. Since, as emphasized in Sec. 2.4, the gravitational 
path integral only depends in a non-trivial way on the dimensionless combination G y/X~o (see 
for example Eqs. (2.30) and (2.31) and related discussion), Newton's constant itself G can be 
decomposed uniquely into non-running and running parts, in the following way 



1 („ a^\2 1 



G = G\- 'y GVX °j e 



(GVAo)Gu 2 ) " (2.51) 



2 



where the running of the second term can be directly deduced from either Eqs. (2.22) or (2.21) 
(both only depend on the combination G\f\~o), or Eq. (2.10), which is related to the previous two 
by the scaling assumption for the free energy F. 

In conclusion, the modified Einstein equations, incorporating the proposed quantum running of 
G, read 

iV - \9v»R + *9nu = 8irG{a)T^ (2.52) 

with A ~ |s, an d only G(y?) on the r.h.s. scale-dependent. The precise meaning of G(0) is given 
in section 4. . 

2.6 Non-perturbative gravitational vacuum condensate 

In this section we will point out the deep relationship (well understood in strongly coupled non- 
abelian gauge theories) between the non-perturbative scale £ appearing in Eqs. (2.1) and (2.2), and 
the non-perturbative vacuum condensate of Eqs. (2.24) and (2.48), which is a measure of curvature. 
The principal, and in our opinion inescapable, conclusion of the results of Eqs. (2.24) and (2.47) is 
that the scale £ appearing in Eqs. (2.1) and (2.2) is related to curvature, and must be macroscopic 
for the theory to be consistent. How can quantum effects propagate to such large distances and 
give such drastic modifications to gravity? The answer to this paradoxical question presumably lies 
in the fact that gravitation is carried by a massless particle whose interactions cannot be screened, 
on any length scale. 

It is worth pointing out here that the gravitational vacuum condensate, which only exists in 
the strong coupling phase G > G c , and which is proportional to the curvature, is genuinely non- 
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perturbative. If one uses a shorthand notation 1Z for it, then one can summarize the result of 
Eq. (2.48) as 

n ~ (l(T 30 ey) 2 ~ c 2 (2-53) 

where the condensate is, according to Eqs. (2.10) and (2.26) (relating £ to \G—G C \), and Eqs. (2.47) 
and (2.48) (relating 1Z to £), non-analytic at G = G c , 

ft ~ | C - G c | 2 ^ (2.54) 

The non-perturbative curvature scale £ then corresponds to a non-vanishing graviton vacuum con- 
densate of order ~ 10~ 30 eV, extraordinarily tiny compared to the QCD color condensate 
(Ajgg ~ 220 MeF) and the electro-weak Higgs condensate (v ~ 250 GeV). But as previously em- 
phasized, the quantum gravity theory cannot provide, at least in its present framework, a value for 
the non-perturbative curvature scale £, which ultimately can only be fixed by phenomenological 
input, either by Eq. (2.1) or, equivalently, by Eq. (2.48). The main message here is that the scale 
in those two equations is one and the same. 

Pursuing the analogy with strongly coupled Yang-Mills theories and QCD, we note that there 
the non-perturbative gluon vacuum condensate depends in a nontrivial way on the corresponding 
confinement scale parameter [40], 

a s < ■ F^ v > ~ (250M e y) 4 ~ £~ 4 (2.55) 

with £qq D ~ Ajjg. The above condensate is not the only one that appears in QCD, another 
important non-vanishing vacuum condensate being the fermionic one [42] 

(a 5 ) 4//3 ° < i>i) > *t - (230M e y) 3 ~ r 3 (2-56) 
2.7 Quantum gravity near two dimensions 

The result of Eqs. (2.1) and Eq. (2.9) are almost identical in overall structure to what one finds 
in 2 + e dimensions, if one allows for a different value of exponent v as one transitions from two 
dimensions to the physical case of four dimensions. In this section we will explore their relationship, 
and the lessons one learns from similar field theory models, such as the non-linear sigma model above 
two dimensions. But the second major ingredient which non-perturbative lattice studies provide 
[20, 26], besides the existence of a phase transition between two geometrically rather distinct phases, 
is that the weakly coupled small G phase is pathological, in the sense that the theory becomes 
unstable, with the four-dimensional lattice collapsing into a tree-like two-dimensional structure for 
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G < G c . Indeed the lattice theory close to the transition at G c is "on the verge" of becoming two- 
dimensional [26], but only to the extent that the effects of higher derivative terms and conformal 
anomaly contributions can be ignored at short distances 6 . 

To one loop the d = 2 + e result for the gravitational beta function was computed some time 
ago, and reads 

/3(G) = (d - 2) G - fa G 2 - . . . , (2.58) 

It exhibits the celebrated non-trivial ultraviolet fixed point of 2 + e quantum gravity at G c = 
(d — 2)/(3o (one finds /3q > for pure gravity). Furthermore, the physics of it bears some striking 
similarity to the non-linear sigma model, to be discussed below. The latter is also perturbatively 
non-renormalizable above two dimensions, but can be constructed by a suitable double expansion 
in the coupling g and e = d — 2 [48] . 

In gravity the corresponding running of G is obtained by integrating Eq. (2.58) 

0(* 2 ) - i ± { J^ r -«n < 2 - 59 > 

The choice of + or — sign is determined from whether one is to the left (+), or to right (-) of the 
fixed point, in which case the effective G(k 2 ) decreases or, respectively, increases as one flows away 
from the ultraviolet fixed point at G c . Physically they represent a screening and an anti-screening 
solution. It is noteworthy that the invariant mass scale m arises as an arbitrary integration constant 
of the renormalization group equations. While in the continuum both phases, and therefore both 
signs, seem acceptable (giving rise to both a "Coulomb" or "spin wave" phase, and a strong coupling 
phase), the euclidean lattice rules out the small G < G c phase as pathological, in the sense that 
the lattice collapses into a two-dimensional branched polymer [20, 13] 7 . 

Thus the smooth phase with G > G c emerges as the only physically acceptable phase in d = 3 
and d = 4 [24, 13]. These arguments suggest therefore that in the renormalization group solution 
Eq. (2.59) only the - sign is meaningful and physical, corresponding to an infrared growth of the 

6 One way of determining coarse aspects of the underlying geometry is to compute the effective dimension in the 
scaling regime, for example by considering how the number of points within a thin shell of geodesic distance between 
t and t + A scales with the geodesic distance itself [26]. For distances a few multiples of the average lattice spacing 
one finds 

N(t) ~ r d \ (2.57) 

T — ► OO 

with d v — 3.1(1) for G > G c (the smooth phase) and d v — 1-6(2) for G < G c (the rough phase). One concludes that 
in the rough phase the lattice tends to collapse into a degenerate tree-like configuration, whereas in the smooth phase 
the effective dimension of space-time is consistent with four. Higher derivative terms tend to affect these results at 
very short distances, where they tend to make the geometry smoother [20, 47]. 

7 The collapse stops at d — 2 because the gravitational action becomes a topological invariant. 
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coupling for G > G c , the gravitational anti-screening solution given by 

2\ (d-2)/2 



G{k 2 ) ~ G. 



c 




+ ... 



(2.60) 



The above expression has in fact exactly the same structure as the lattice result of Eq. (2.1). 
Even physically the result makes sense, as one would expect that gravity cannot be screened (as 
would happen for the "+" sign choice). Repeating some of the general arguments presented in the 
previous sections, and replacing k 2 — ► A 2 in the above equation, where A is the ultraviolet cutoff, 
one can solve, for G > G c , for the genuinely non-perturbative, dynamically generated mass scale 
m in terms of the coupling at the cutoff (lattice) scale, as in Eq. (2.16). It should be noted here 
that Eq. (2.16) is essentially the same as Eq. (2.10), and that Eq. (2.59) is essentially the same as 
Eq. (2.1). 

The derivative of the beta function at the fixed point gives the exponent v. From Eq. (2.58) 
one has therefore close to two dimensions [44] 

P'{G C ) = -l/u = (2 - d) (2.61) 

Recently the above results have been extended to two loops, giving close to two dimensions 8 

2 20 

(3(G) = {d - 2) G - - (25 - n s ) G 2 — (25 - n s )G 3 + . . . , (2.62) 

for n s massless real scalar fields minimally coupled to gravity. After solving the equation (3(G C ) = 
to determine the location of the ultraviolet fixed point, one finds 



_3 ( t} _o\ 

2 (25 - n s ) ^ ~ >~ 2(25 - n s 



°c = 0/oc — (d-2)- — : (d-2) 2 + ... 



v 



-l nUn \ (a o \ i i a o^2 



0(G C ) = (d - 2) + (d - 2) 2 + . . . (2.63) 

Tig 

which gives, for pure gravity without matter [n s = 0) in four dimensions, to lowest order v~ x = 2, 
and v~ l 4.4 at the next order. Unfortunately in general the convergence properties of the 2 + e 
expansion for some better understood field theories, such as the nondinear sigma model, are not 
too encouraging, at least when compared for example to well-established results obtained by other 
means directly in d = 3 [51, 52]. This somewhat undesirable state of affairs is usually ascribed to 

8 For a while there was considerable uncertainty, due in part to the kinematic singularities which arise in gravity 
close to two dimensions, about the value of the graviton contribution to /3o, which was quoted originally in [43] as 
38/3, in [44] as 2/3, and more recently in [45, 46] as 50/3. As discussed in [8], the original expectation was that 
the graviton contribution should be d(d — 3)/2 = —1 times the scalar contribution close to d = 2. Direct numerical 
estimates in d — 3 give v~ x ~ 1.67 [24] and are therefore in much better agreement with the larger, more recent value 
for /?o. 
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the suspected existence of renormalon-type singularities ~ e~ c ^ G close to two dimensions, which 
could possibly arise in gravity as well. At the quantitative level, the results of the 2 + e expansion 
for gravity remain therefore somewhat limited, and obtaining the three- or four-loop corrections 
could represents a daunting task. But one should not overlook the fact that there are quantum 
field theory models which have in some qualitative respects astonishingly similar behavior in 2 + e 
dimensions, and are at the same time very well understood. Some analogies might therefore remain 
helpful. 

The non-linear O(N) sigma model is one example, studied extensively in the context of the 
2 + e expansion, and solved exactly in the large N limit [38]. The model is not perturbatively 
renormalizable above two dimensions. Yet in both approaches it exhibits a non-trivial ultraviolet 
fixed point at g c (a phase transition in statistical mechanics language), separating a weak coupling 
massless ordered phase from a massive, strong coupling phase. Therefore the correct continuum 
limit has to be taken in the vicinity of the non-trivial ultraviolet fixed point. Perhaps one of the 
most striking aspects of the non-linear sigma model above two dimensions is that all particles are 
massless in perturbation theory, yet they all become massive in the strong coupling phase g > g c , 
with masses proportional to a non-perturbative scale m [49]. 

The second example of perturbatively non-renormalizable theory is the chirally-invariant self- 
coupled fermion above two dimensions [38]. It too exhibits a non-trivial ultraviolet fixed point 
above two dimensions, and can be studied perturbatively via a double 2 + e expansion. And it too 
can be solved exactly in the large N limit. With either method, one can show that the model is 
characterized by two phases, a weak coupling phase where the fermions stay massless and chiral 
symmetry is unbroken, and a strong coupling phase in which chiral symmetry is spontaneously 
broken, a fermion condensate arises, and a mass scale is generated non-perturbatively. 

2.8 Other determinations of the exponent v 

Let us conclude this discussion by mentioning the remaining methods which have been used to 
estimate v. 

Recently some approximate renormalization group results have been obtained in the continuum 
based on an Einstein-Hilbert action truncation. In the limit of vanishing bare cosmological constant 
the result u^ 1 = 2d(d — 2)/(d + 2) = 2.667 was given in d = 4 [53]. In the cited work the sensitivity 
of the numerical answer for the exponent to the choice of gauge fixing term and to the specific shape 
of the momentum cutoff were systematically investigated as well. More details about the procedure 
can be found in the quoted references. The more recent results extend earlier calculations for the 
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exponent v done by similar operator truncation methods, described in detail in [55] and references 
therein. A quantitative comparison of these continuum results with the lattice answer for is^ 1 , 
and with the 2 + e estimates, can be found in Figs. 1 and 2. The overall the agreement seems 
reasonable, given the uncertainties inherent in the various estimates. There is also some connection 
between the just mentioned lattice and continuum renormalization group results and the ideas in 
[56], but in this last reference fractional exponents, such as the ones in Eq. (2.2) and (2.59), were 
not considered. 

Both the lattice and the continuum renormalization group analysis can be extended to dimen- 
sions greater than four. In [14] it was suggested, based on a simple geometric argument, that 
v = l/(d — 1) for large d, and in the large-d limit one can prove that v = [15]. Moreover, for the 
lattice theory in finite dimensions one finds no phase transition in d = 2 [23], v ~ 0.60 in d = 3 
[24] and v ~ 0.33 in d = 4 [13, 14], which leads to the more or less constant sequence {d — 2)u = 
1, 0.60 and 0.66 in the three cases respectively, which would suggest v ~ 1/d for large d. On the 
other hand, in the same large d limit, in [53] the value v ~ l/2d was obtained, again performing a 
renormalization group analysis truncated to the Einstein-Hilbert action with a cosmological term 
in the continuum. One cannot help noticing some reasonable agreement between the estimates for 
v coming from the lattice, and the corresponding results for v using the truncated renormalization 
group expansion in the continuum, as well as some degree of compatibility with the 2 + e expansion 
close to d = 2. For a more quantitative comparison, see again Figs. 1 and 2. 

One possible advantage of the continuum renormalization group calculations is that, since they 
can be performed to a great extent analytically, they allow greater flexibility in exploring various 
scenarios, such as additional invariant operators in the action, measure contributions, or varying 
dimensions. At the same time the absence of a reliable estimate for the errors involved in the 
truncation leaves the method with uncertainties, which might be hard to quantify until an improved 
calculation is performed with an extended operator basis. 
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Fig. 1. Universal gravitational exponent \jv of Eqs. (2.1) and (2.2) as a function of the dimension (the 
abscissa here is z = (d — 2)/{d— 1), where d is the space-time dimension, which maps d = 2 to z — and 
d = oo to z = 1). The larger circles at g? = 3 and d = 4 are the lattice gravity results of [24, 13], interpolated 
(continuous curve) using the exact lattice results 1/v — in d — 2 [23], and v — at d = oo [15]. The smaller 
black dots (connected by the dashed curve) are the recent results of [53, 54], obtained from a continuum 
renormalization group study around the non-trivial fixed point of the gravitational action in d dimensions 
in the limit A —> 0. The two curves close to the origin are the 2 + e expansion for 1/v to one loop (lower 
curve) and two loops (upper curve) [46] . 
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Fig. 2. Critical point k c = l/(8irG c ) in units of the cutoff (as it appears for ex. in Eqs. (2.10) and (2.22)) 
as a function of dimension. The abscissa is the space-time dimension d. The large circles at d = 3 and 
d = 4 are the Regge lattice results of [24, 13], interpolated (dashed curve) using the additional lattice results 
l/k c = in d = 2 [23]. The lower continuous curve is the large-d lattice result of Eq. (2.27). The two 
upper continuous curves are the two-loop 2 + e expansion result in the continuum of Eq. (2.63) [46], whose 
prediction becomes uncertain as d approaches the value three. 

2.9 The running of a(/j, 2 ) in QED and QCD 

QED and QCD provide two invaluable illustrative cases where the running of the gauge coupling 
with energy is not only theoretically well understood, but also verified experimentally. This section 
is intended to provide analogies and distinctions between the two theories, and later gives justifica- 
tion, based on the structure of infrared divergences in QCD, for the transition between the running 
of the gravitational coupling as given in Eq. (2.1), and its infrared regulated version of Eq. (2.2). 
Most of the results found in this section are well known, but the purpose here is provide a contrast 
(and in some instances, a close relationship) with the gravitational case. 

In QED the non-relativistic static Coulomb potential is affected by the vacuum polarization 
contribution due to electrons (and positrons) of mass m. To lowest order in the fine structure 
constant, the contribution is from a single Feynman diagram involving a fermion loop. One finds 
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for the vacuum polarization contribution u;^(k 2 ) at small k 2 the well known result [57] 

(2.64) 



e 2 e 2 e 2 



a k 2 9 

^ + ° « 

15 7T m z 



k 2 k 2 [l+c^(k 2 )] k 2 

which, for a Coulomb potential with a charge centered at the origin of strength — Ze leads to 
well-known Uehling [58, 59] 5- function correction 

v{r) = ( 1 _^A)z^ = z^_^z^! J( 3, (x) (M6) 

V 157rm 2 / Airr Airr I5ir m 2 

It is not necessary though to resort to the small-k 2 approximation, and in general a static charge 
of strength e at the origin will give rise to a modified potential 

Q(r) (2.66) 



47rr 4irr 
with 



a 1 

Q(r) = 1 + — In —5-5 + ... mr«l (2.67) 



3 7r m 2 r 2 
for small r, and 

3(0 = 1 + , ° W 2 e " 2mr + ••• ( 2 - 68 ) 
4 ^/7T (mr)' i / z 

for large r. Here the normalization is such that the potential at infinity has Q(oo) = 1 9 . The reason 
we have belabored on this example is to show that the screening vacuum polarization contribution 
would have dramatic effects in QED if for some reason the particle running through the fermion 
loop diagram had a much smaller (or even close to zero) mass. 

There are two interesting aspects of the (one- loop) result of Eqs. (2.67) and (2.68). The first one 
is that the exponentially small size of the correction at large r is linked with the fact that the electron 
mass m e is not too small: the range of the correction term is £ = 2h/mc = 0.78 x 10 _10 cm, but 
would have been much larger if the electron mass had been a lot smaller. In fact in the limit of zero 
electron mass the correction progressively extends out all the way to infinity. Ultimately the fact 
that the Uehling correction is not important in atomic physics is largely due to the fact that its value 
is already quite small on atomic scales, comparable to the Bohr radius ao = h/me 2 = 0.53 x 10~ 8 . 

The second interesting aspect is that the correction is divergent at small r due to the log term, 
giving rise to a large unscreened charge close to the origin; and of course the smaller the electron 
mass m, the larger the correction. Also, it should be noted that the observed laboratory value for 

the effective electromagnetic charge is normalized so that Q(r) is Q(r = oo) = 1, whereas in gravity 

9 The running of the fine structure constant has recently been verified experimentally at LEP, the scale dependence 
of the vacuum polarization effects gives a fine structure constant changing from a(0) ~ 1/137.036 at atomic distances 
to about a(mz ) ~ 1/128.978 at energies comparable to the Z boson mass, in good agreement with the theoretical 
prediction [60]. 
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the laboratory value corresponds to the "short distance" limit, G(r = 0), due to the very large size 
off in Eq. (2.2). 

In QCD (and related Yang-Mills theories) radiative corrections are also known to alter signif- 
icantly the behavior of the static potential at short distances. The changes in the potential are 
best expressed in terms of the running strong coupling constant as(fi), whose scale dependence is 
determined by the celebrated beta function of SU(3) QCD with nf light fermion flavors [61] 



»-ej = = - j^as ~ ^3«5 - ••• (2-69) 



with Pq = 11 - |n/, [3i = 51 - fn f , and (3 2 = 2857 - ^n f + ^-nj. The solution of the 
renormalization group equation Eq. (2.69) then gives for the running of as(n) 



4 7T 

as(fi) 



_ 2A ln[ln^/Af^] 



2 _ 

2.70) 



The non-perturbative scale A-^ appears as an integration constant of the renormalization group 
equations, and is therefore - by construction - scale independent. The physical value of A^jg cannot 
be fixed from perturbation theory alone, and needs to be determined by experiment 10 . 

In QCD one determines experimentally from Eq. (2.70) (for example, from the size of scaling 
violations in deep inelastic scattering) A-p^ ~ 220MeV, which not surprisingly is close to a typical 
hadronic scale. But from a purely theoretical standpoint one could envision other scenarios, where 
this scale would take on a completely different value. The point is that nothing in QCD itself 
determines the absolute magnitude of this scale: only its size in relation to other physical observables 
such as hadron masses is theoretically calculable within QCD. 

In principle one can solve for Aj^g in terms of the coupling at any scale, and in particular at 
the cutoff scale A, obtaining 

A M = Aexp^-y = A (^fp) e-^m [1 + 0(a s (A))} (2.71) 

It should be clear by now that this expression is the analog of Eqs. (2.10) and (2.16) in the 
gravitational case. In lattice QCD this is usually taken as the definition of the running strong 
coupling constant as (/x). It then leads to an effective potential between quarks and anti-quarks of 
the form 

V(k') = -12^1 (2.72) 



10 Alternatively, the integration constant in Eq. (2.69) can be fixed by specifying the value of as at some specific 
energy scale, usually chosen to be the Z° mass. 
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and the leading logarithmic correction makes the potential appear softer close to the origin, V(r) ~ 
l/(r lnr). 

When the QCD result is contrasted with the QED answer of Eqs. (2.64) and (2.65) it appears 
that the infrared small k 2 singularity in Eqs. (2.72) is quite serious. An analogous conclusion is 
reached when examining Eqs. (2.70) : the coupling strength as(k 2 ) diverges in the infrared due 
to the singularity at k 2 = 0. In order to avoid a meaningless divergent answer, the uncontrolled 
growth in as(k 2 ) needs therefore to be regulated (or self-regulated) by the dynamically generated 
QCD infrared cutoff A-^r [64] (as indeed happens in simpler theories, such as the non-linear sigma 
model, which can be solved exactly in the large N limit, where this mechanism can therefore be 
explicitly demonstrated [38]). To lowest order in the coupling this implies 



2 I 47T i 47T 



Vfr) = t s- 7 s- (2.73) 

3 0ok 2 ln ( k2/A y 3 Aka ln ( (k2 + A 2_ } /A y 



The resulting potential can then be evaluated for small k 2 
giving, after Fourier transforming back to real space, for large r 



y <*> *,«L« -zww* - (2 ' 74) 



V{r) -~=o "t^W' (-8V'') + ^ = + " (2 ' 75) 

The desired linear potential, here with string tension a, is then indeed recovered at large distances. 
In fact the interpolating potential of Eq. (2.73) is remarkably successful in describing non-relativistic 
QCD bound states, as discussed in [64, 65], incorporating correctly to some extent the leading short- 
and large-distance QCD corrections 11 . 

What is rather remarkable in this context is that the removal of a serious infrared divergence in 
as(k 2 ), by the replacement of Eq. (2.73), has in fact caused an even stronger infrared divergence in 
y(k 2 ) itself, which leads though, in this instance, to precisely the desired result of Eq. (2.75), namely 
a confining potential at large distances, as expected on the basis of non-perturbative arguments 
[35]. From as (k 2 ) given in Eq. (2.73) one can calculate the corresponding beta function 



P(*s) = -^a 2 s [l-e eo°s) + 0(a s ) (2.76) 

which now exhibits a non-analytic (renormalon) contribution at as = 0, not detectable to any 
order in perturbation theory [62, 63]. 

11 A similar infrared regularization was done in Eq. (2.2), (m 2 /k 2 ) 1 / 2l/ = exp [— ln(k 2 /m 2 )] — » 
cxp{^lri[(k 2 +m 2 )/m 2 ]} 
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The very fruitful analogy between strongly coupled non-abelian gauge theories and strongly 
coupled gravity can be pushed further [15], by developing other aspects of the correspondence 
which naturally lend themselves to such a comparison. In gravity the analog to the Wilson loop 
W(r) of non-abelian gauge theories 



W(T) ~ Tr V exp 



An dX^ 

c 



(2.77) 



exists as well (see Eq. (2.45)), defined in the gravitational case as the path-ordered exponential of 
the affine connection T^ u around a closed planar loop. 

Furthermore, it is known that in QCD, due to the non-trivial strong coupling dynamics, there 
arise several non-perturbative condensates. Thus the gluon condensate is related to the confinement 
scale via [40] 

a s < iV ■ F^ v > ~ (250M e y) 4 ~ £ -4 (2.78) 

with ~ A^. Another important non-vanishing non-perturbative condensate in QCD is the 

fermionic one [41, 42] 

(a 5 ) 4//3 ° < > ^ -(230MeF) 3 ~ C 3 (2-79) 

and represents a key quantity in discussing the size of chiral symmetry breaking, and the light 
quark masses. 



3 Poisson's Equation with Running G 



Given the running of G from either Eq. (2.2), or Eq. (2.1) in the large k limit, the next step is 
naturally a solution of Poisson's equation with a point source at the origin, in order to determine 
the structure of the quantum corrections to the gravitational potential in real space. The more 
complex solution of the fully relativistic problem will then be addressed in the following sections. 
In the limit of weak fields the relativistic field equations 

-\g„vR + \g^ = SnGT^ (3.1) 

give for the cf> field (with goo {%) — — (1 + 2<p(x) ) 

(A - \)<f>[x) = 4ttG P (x) - A (3.2) 

which would suggest that the scaled cosmological constant A acts like a mass term m = \/\. For a 
point source at the origin, the first term on the r.h.s is just AttMG 6^ (x). The solution for <j)(r) 
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can then be obtained simply by Fourier transforming back to real space Eq. (2.6), and, up to an 
additive constant, one has 

</>(r) = -MG— - aomMG ( mr )*(£-i)# ( mr ) (3.3) 

where K n {x) is the modified Bessel function of the second type. The behavior of <j>[r) would then 
be Yukawa-like 0(r) ~ const. e~ mr /r and thus rapidly decreasing for large r. 

But the reason why both of the above results are in fact incorrect (assuming of course the validity 
of general coordinate invariance at very large distances r>l/ \/A) is that the exact solution to the 
field equations in the static isotropic case with a A term gives 

2MG A 2 

-5oo = B(r) = 1 -r 2 (3.4) 

V o 

showing that the A term definitely does not act like a mass term in this context. 

Therefore the zeroth order contribution to the potential should be taken to be proportional to 
47r/(k 2 + /i 2 ) with fi — > 0, as already indicated in fact in Eq. (2.8). Also, proper care has to be 
exercised in providing an appropriate infrared regulated version of G(k 2 ), and therefore V(k 2 ), 
which from Eq. (2.8) reads 



47T 



(k 2 + fi 2 ) 



1 + a 



m 2 \ 2v 



k 2 + m 



2 



(3.5) 



and where the limit p — > is intended to be taken at the end of the calculation. 

There are in principle two equivalent ways to compute the potential 4>(r), either by inverse 
Fourier transform of the above expression, or by solving Poisson's equation A(p = 4np with p(r) 
given by the inverse Fourier transform of the correction to G(k 2 ), as given later in Eq. (3.17). Here 
we will first use the first, direct method. 

3.1 Large r limit 

The zeroth order term gives the standard Newtonian —MG/r term, while the correction in general 
is given by a rather complicated hyper geometric function. But for the special case v = 1/2 one has 
for the Fourier transform of the correction to <j){r) 

i/„ 4vr 1 2 e-^ r -e- mr 2 1 - e~ mr , , 

a m ' -5— x -> a m — — ^— ~ a m ^ (3.6) 

k 2 + p 2 (k 2 + m 2 ) 277 r(m 2 - p 2 ) m^o m 2 r 

giving for the complete quantum-corrected potential 

<K0 = [l + G (l-e- mr )] (3.7) 
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For this special case the running of G(r) is particularly transparent, 

G(r) = Goo ( 1 - -^—e~ mr ) (3.8) 

with Gqq = (1 + 0,0) G and G = G(0). G therefore increases slowly from its value G at small r to 
the larger value (1 + ao) G at infinity. Figure 3. provides a schematic illustration of the behavior 
of G as a function of r. 

3 „ , , , , , , , , , , , , , , , , , , , „ 

2.5 : : 



u 

"1.5- 




0.5- 



5 10 15 20 

r 



Fig. 3. Schematic scale dependence of the gravitational coupling G(r) from Eq. (3.8), here for v = 1/2. 
The gravitational coupling rises initially like a power of r, and then approaches the asymptotic value Goo = 
(1 + ao)G for large r. The behavior for other values of v > 1/3 is similar. 

Returning to the general v case, one can expand for small k to get the correct large r behavior, 
1 1111 



1 

(k 2 + /i 2 )(k 2 + m 2 )277 ~~ m \ (k 2 + /i 2 ) 2um 2 ( k 2 + m 2)^ 

After Fourier transform, one obtains the previous answer for v = 1/2, whereas for v = 1/3 one 
finds 



(3.9) 



MG 



1 + ao 1 



2>mr 



7T 



Ko(mr) 



and for general v 



MG 



l + oo 1 



2 2 ( 3 J mr , , _ i cq i \ 



2v 



(3.10) 



(3.11) 
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Using the asymptotic expansion of the modified Bessel function K n {x) for large arguments, K n (z 
\A"/2 z~ l l 2 e~ z (1 + 0(l/z)), one finally obtains in the large r limit 

M G 



4>{r) 

with Q = 1/(1/ 2k r(^)). 



r — ► oo 



1 + a (l - q(mr)277- 1 e- mr ) 



(3.12) 



(-MG) a m 2 



3vr 



ln(- 



+ 0(r 3 



(3.13) 



3.2 Small r limit 

In the small r limit one finds instead, using again Fourier transforms, for the correction for v = 1/3 

•■- ,mr 5 

In the general case the complete leading correction to the potential <p(r) for small r (and v > 1/3) 
has the structure (— const. )(—MG)ao m~ r~~ l . Note that the quantum correction always vanishes 
at short distances r — ► 0, as expected from the original result of Eqs. (2.1) or (2.2) for k 2 — ► oo 12 . 

The same result can be obtained via a different, but equivalent, procedure, in which one solves 
directly the radial Poisson equation for 4>{r). First, for a point source at the origin, x , 

with 



tf(3) (x) = J_ 2 



«5(r) 



4-7T 



(3.14) 



one sets A(f>(r) — > r -^/dr^r^r)] in radial coordinates. In the ao / case one then needs to 
solve Acp = 4np, or in the radial coordinate for r > 



1 d 



di 



4:TrGp m (r) 



(3.15) 



m 

k 2 + m 2 



2v 



r 2 dr V 

with the source term p m determined from the inverse Fourier transform of the correction term in 
Eq. (2.2), namely 

One finds 

Pm(r) = -^-c„a Mm 3 (mr)"^ 3 "^ Ki,„_iJmr) 

o7T 2V° i/> 

with 

c u = 



v^r(^) 

The vacuum polarization density p m has the property 



(3.16) 
(3.17) 
(3.18) 



47T / r 2 dr p m (r) = a 
Jo 



M 



(3.19) 



12 At very short distances r ~ ip other quantum corrections come into play, which are not properly encoded in 
Eq. (2.1), which after all is supposed to describe the universal running in the scaling region lp <JC r <C £. Furthermore, 
higher derivative terms could also have important effects at very short distances. 
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where the standard integral f£° dx x 2 ~ n K n (x) = 2 _n y / 7rr ^| — nj has been used. Note that the 
vacuum polarization distribution is singular close to r = 0, just as in QED, Eq. (2.67). 

The r — > result for (j>{r) (discussed in the following, as an example, for v = 1/3) can then be 
obtained by solving the radial equation for (f>(r), 



_ r 0( r ) = K (mr) 3.20 

r dr z it 



where the (modified) Bessel function is expanded out to lowest order in r, Ko(m r) = —7— In (^j 1 ) + 
0(m 2 r 2 ), giving 



2 r 



MG 3 r 

cj)( r ) = + a G MGm 6 — 

r 37r 



, ,mr. 5 
-ln(— )-7+e 



+ 0(r 3 ) (3.21) 



where the two integration constants are matched to the large r solution of Eq. (3.11). Note again 
that the vacuum polarization density p m (r) has the expected normalization property 

f°° o , a M m 3 „ . , 2a Mm 3 it 
4vr / r 2 dr " g mr = — — = a M 3.22 







so that the total enclosed additional "charge" is indeed just aoM, and = Go(l + Oq) (see for 
comparison also Eq. (3.11)). Using then the same method for general v > |, one finds for small r 
(using the expansion of the modified Bessel function K n (x) for small arguments as given later in 
Eq. (5.24)) 

Pm{r) ~ |SeC / , x Q o M mv r^ 3 = A rH» (3.23) 
4vrr(l-l) 

and from it the general result 

4>(r) ~ + ao MGc s mv r^ 1 + ... (3.24) 

r— > r 

with c a = i/|sec(|;)|/r(i). 

4 Relativistic Field Equations with Running G 



Solutions to Poisson's equation with a running G provide some insights into the structure 
of the quantum corrections, but a complete analysis requires the study of the full relativistic 
field equations, which will be discussed next in this section. A set of relativistic field equations 
incorporating the running of G is obtained by doing the replacement [12] 

G(k 2 ) - <?(□) (4.1) 
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with the d'Alembertian □ intended to correctly represent invariant distances, and incorporating 
the running of G as expressed in either Eqs. (2.1) or (2.2), 



G -► G(u) = G 



1 + 00 + - 



(4.2) 



For the use of □ to express the running of couplings in gauge theories the reader is referred to the 
references in [66]. Here the □ operator is defined through the appropriate combination of covariant 
derivatives 

□ = gTVuVv (4.3) 

and whose explicit form depends on the specific tensor nature of the object it is acting on, as in 
the case of the energy-momentum tensor 

D T aP - jS ... = sT V„ (V„ T°% ) (4.4) 
Thus on scalar functions one obtains the fairly simple result 

nS(x) = -^d^g^^d u S(x) (4.5) 
V 9 

whereas on second rank tensors one has the significantly more complicated expression □ T a/ g = 
gV v V M (V^r Q/ g). In general the invariant operator appearing in the above expression, namely 

= a {—) (4.6) 

or its infrared regulated version 

has to be suitably defined by analytic continuation from positive integer powers; the latter can 
be often be done by computing O n for positive integer n, and then analytically continuing to 
n — > —l/2u. In the following the above analytic continuation from positive integer n will always 
be understood. Usually it is easier to work with the expression in Eq. (4.6), and then later amend 
the final result to include the infrared regulator, if needed. 

One is therefore lead to consider the effective field equations of Eq. (2.52), namely 

iV - \9^R + A<V = 8vrG (l + A(U)) (4.8) 

with -A(D) given by Eq. (4.7) and A ~ l/£ 2 , as well as the trace equation 

R - 4A = -8ttG (1 + A(n)) T (4.9) 

34 



Being manifestly covariant, these expressions at least satisfy some of the requirements for a set of 
consistent field equations incorporating the running of G, and can then be easily re-cast in a form 
similar to the classical field equations 

iV -\g^R + \g^ = SnGf^ (4.10) 

with T^v = (1 + A(U)) T^ v defined as an effective, gravitationally dressed, energy-momentum ten- 
sor. Just like the ordinary Einstein gravity case, in general might not be covariantly conserved 
a priori, V M T^ v / 0, but ultimately the consistency of the effective field equations demands that 
it be exactly conserved in consideration of the Bianchi identity satisfied by the Einstein tensor [12]. 
The ensuing new covariant conservation law 

V%„ = V [(l + A(U)) T„ v ] = (4.11) 

can be then be viewed as a constraint on T^ v (or T^ v ) which, for example, in the specific case of 
a perfect fluid, implies a definite relationship between the density p(t), the pressure p(t) and the 
metric components [12]. 

From now on, we will set the cosmological constant A = 0, and its contribution can then be 
added at a later stage. As long as one is interested in static isotropic solutions, one takes for the 
metric the most general form 

ds 2 = -B{r)dt 2 + A{r)dr 2 + r 2 (d6 2 + sin 2 9 dip 2 ) (4.12) 

and for the energy momentum tensor the perfect fluid form 

T^ u = diag [B(r) p(r), A(r)p(r), r 2 p(r), r 2 sin 2 6p(r) ] (4.13) 

with trace T = 3p — p. The trace equation then only involves the (simpler) scalar d'Alembertian, 
acting on the trace of the energy-momentum tensor 

To the order one is working here, the above effective field equations should be equivalent to 

I + \ {a) (R^-\%vR) = (1 - A(U) + = SirGT^ (4.14) 

where the running of G has been moved over to the "gravitational" side, and later expanded out, 
assuming the correction to be small. For the vacuum solutions, the r.h.s. is zero for r / 0, and one 
can re-write the last equation simply as 

1 



8vrG(l + A{n)) 



(R^-^g^R) = (4.15) 
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4.1 Dirac delta function source 

A mass point source is most suitably described by a Dirac delta function. The delta function at 
the origin can be represented for example as 

6(r) = lim ——^ — ~r (4.16) 

w £ -o 7T (r 2 + e 2 ) ; 

As an example, the derivative operator (d/dr) n acting on the delta function would give, for small 
e, 

(-l) n r(n + 2)vr- 1 er n (r 2 + e 2 )""" 1 (4.17) 
which can be formally analytically continued to a fractional value n = — 1/v 

(-l)-V- r (^2 - I) .-^r 1 ^ ( r 2 + e 2 ) (4.18) 
and then re-expressed as an overall factor in front of the original 5 function 

(-l)- 1 /^- 1 ^ (r 2 + e 2 )^ T 5(r) (-lj-VV/T^-^j <5(r) (4.19) 

The procedure achieves in this case the desired result, namely the multiplication of the original 
delta function by a power r x l v . Note also that if one just takes the limit e — > at fixed r one always 
gets zero for any r > 0, so close to the origin r has to be sent to zero as well to get a non-trivial 
result. 

The above considerations suggests that one should be able to write for the point source at the 
origin 

T^ir) = diag[B(r)p(r), 0, 0, 0] (4.20) 

with (in spherical coordinates) 

p = M^(x) -► (4.21) 

47r r 2 

and the delta function defined as a suitable limit of a smooth function, and with vanishing pressure 
p(r) = 0. 

Next we will consider as a warmup the trace equation 

R = -8irG(D)T = -8ttG (1 + -A(D)) T = +8ttG (1 + p (4.22) 

where we have used the fact that the point source at the origin is described just by the density 
term. One then computes the repeated action of the invariant d'Alembertian on T, 

, „ , 16G7rp' AGnA'p' AGnB'p' 8Gtt/ Ia . 

□ (-SvrGT) = D^Cp) = ^ + Ag + — f- (4.23) 
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or, using the explicit form for p(r), 

□ (_ 8 , T) _ - 8CM ' ( - 6 f- i ;^: riB ' ) + Of) (424) 

7r(r 2 + e 2 ) 3 A 2 B 

In view of the rapidly escalating complexity of the problem it seems sensible to expand around the 
Schwarzschild solution, and set therefore 

A(r)-' = 1 - ™° + °M (4.25) 



and 



2MG 9{r) , A , 

B(r) = 1 + — 4.26 

r r 



where the correction to the standard solution are parametrized here by the two functions a(r) 

and 9{r), both assumed to be "small", i.e. proportional to ao as in Eqs. (4.6) or (4.7), with ao 

considered a small parameter. Then for the scalar curvature, to lowest order in a and 9, one has 

GM (a — 9) — (2G M — r) (G M 9' + (3GM - 2r) a' + {2GM -r) r9") 

r 2 (r -2MG) 2 ( ' ' 

To simplify the problem even further, we will assume that for 2MG <C r <C £ (the "physical" 

regime) one can set 

a ( r ) = - ao MGc a r a (4.28) 

and 

0(r) = -a MGc e rP (4.29) 

This assumption is in part justified by the form of the non-relativistic correction of Eqs. (3.13). 
Then for a = (3 (the equations seem impossible to satisfy if a and [3 are different) one obtains for 
the scalar curvature 

R = + a (2 c a + (a - 1) eg) a M G r a ~ 3 + +0(a 2 1 ) (4.30) 

A first result can be obtained in the following way. Since in the ordinary Einstein case one has for 
a perfect fluid R = —8nGT = +87rG(p — 3p), and since p m (r) ~ r~~ 3 from Eq. (3.17) in the same 
regime, one concludes that a solution is given by 

a = i (4.31) 

which also seems consistent with the Poisson equation result of Eq. (3.24). 

Next one needs the action of D n on the point source (here hidden in T). To lowest order one 
has for the source term 

8irGp = MG-^-4 ^ ( 4 - 32 ) 
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The d'Alembertian then acts on the source term and gives 



□"(SttGp) - (2n 2 +2)! M(?lr- 2 "- 3 (1) (4.33) 
which can then be analytically continued to n = — resulting in 

(□)-£(8ttGp) r(3 ~ ^ MG^r^' 3 Q (4.34) 

After multiplying the above expression by a , consistency with l.h.s. of the trace equation, 
Eq. (4.30), is achieved to lowest order in ao provided again a = l/u. To zeroth order in ao, 
the correct solution is of course already built into the structure of Eqs. (4.25) and (4.26). Also note 
that setting e = would give nonsensical results, and in particular the effective density would be 
zero for r ^ 0, in disagreement with the result of Eq. (3.17), p m {r) ~ r~~ for small r. 

The next step up would be the consideration of the action of □ on the point source, as it appears 
in the full effective field equations of Eq. (4.8), with again T^ v described by Eq. (4.20). One perhaps 
surprising fact is the generation of an effective pressure term by the action of □, suggesting that 
both terms should arise in the correct description of vacuum polarization effects, 

pB' 2 2Bp' BA'p' gV Bp" 
(OT^) tt - -^^ + — __ + __ + __ 

oB' 2 

(□7» rr = (4.35) 

and (OT F ) M = (OT^,,) = 0. A similar effect, namely the generation of an effective vacuum 
pressure term in the field equations by the action of □, is seen also in the Robertson- Walker metric 
case [12]. 

4.2 Effective trace equation 

To check the overall consistency of the approach, consider the set of effective field equations that 
are obtained when the operator (1 + -A(n)) appearing in Eqs. (4.8) and (4.9) is moved over to the 
gravitational side, as in Eq. (4.15). Since the r.h.s of the field equations then vanishes for r ^ 0, 
one has apparently reduced the problem to one of finding vacuum solutions of a modified, non-local 
field equation. 

Let us first look at the simpler trace equation, valid again for r / 0. If we denote by 5R 
the lowest order variation (that is, of order ao) in the scalar curvature over the ordinary vacuum 
solution R = 0, then one needs to find solutions to 

— 77-r, ^ SR = (4.36) 
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On a generic scalar function F{r) one has the following action of the covariant d'Alembertian □ : 

A'F' B'F' 2F' F" lA „. 

The Ricci scalar is complicated enough, even in this simple case, and equal to 

B' 2 A'B' 2B' 2A' B" 2 2 

O 42 D ~~ Z A D ^ Z A2 ~A~D ~~ Z3, A ~J2 (4.38) 



2AB 2 2A 2 5 r^B rA 2 AB r 2 A r 2 
The action of the covariant d'Alembertian on it produces the rather formidable expression 



5B' 4 3A'B' 3 5B' 3 3A' 2 B' 2 6A'B' 2 5A"B' 2 6B"B' 2 B 
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2A 2 B 4 A 3 B 3 rA 2 B 3 A^B 2 rA 3 B 2 AA 3 B 2 A 2 B 3 r 2 A 2 B 2 

7A' 3 B' 9A' 2 B' A'B' 13A'A"B' AA"B' 23 A' B" B' 9B"B' A^B' 
+ 2A 5 B ~ rA 4 B ~ r 2 A 3 B 4A 4 B + rA 3 B ~ 4A 3 B 2 + rA 2 B 2 + 2A 3 B 

2B' 2B' UA' 3 A A' 2 2B" 2 2 A' 6 A' 2 A" 



+ , .o„ + ^~, o~TT ~t~ ~Tn~^n ~t~ o To o To" ~t~ 



2A 2 B 2 r 3 AB r 3 A 2 B rA 5 r 2 A 4 A 2 B 2 r 3 A 2 r 3 A 3 r 2 A 3 
19A' 2 B" 8A'B" 2A"B" 2A^ 3A'B^ AB^ B^ 4 4 
4A 4 B + 'TAW + A 3 B + TAf + A 3 B ~ TWb ~ A^B + ZT A~ ZT A 2 ( ' 39) 

To lowest order in the functions a and 9, from Eq. (4.27), the scalar curvature is given by 

(GM (a -9)- (2GM - r) (GM9' + (3GM - 2r) a' + (2GM - r) rd")) / (r 2 (r - 2MG) 2 ) (4.40) 

After having the d'Alembertian □ act on this expression, one obtains the still formidable (here 
again to lowest order in a and 9) result 

(9^r 7 + 2<7< 3 V 6 - 8GM9^r 6 - 4a"r 5 + GM9 i3) r 5 - 15GMa^r 5 + 24G 2 mV 4 V 5 + 4a'r 4 
+3GM9"r A + 3lGMa"r A - 6G 2 M 2 9 (3) r 4 + 42G 2 MV 3) r 4 - 32G 3 M 3 9 {A) r A - 12GM9 'r 3 
-2SGMa'r 3 - UG 2 M 2 9"r 3 - 94G 2 M 2 a" r 3 + 12G 3 M 3 9 i - 3) r 3 - 52G 3 M 3 a {3) r 3 + 16G 4 M 4 # (4) r 3 
+ l2GM9r 2 - \2GMor 2 + 48G 2 M 2 0'r 2 + 96G 2 M 2 a'r 2 + 20G 3 M 3 9"r 2 + \32G 3 M 3 a"r 2 
-8G 4 M 4 6> (3) r 2 + 24G 4 MV 3) r 2 - 24G 2 M 2 9r + 2AG 2 M 2 ar - UG 3 M 3 9'r - lQ0G 3 M 3 a'r 
-8G 4 MW - 72G 4 MW + 16G 3 M 3 - 16G 3 M 3 ct + 32G 4 M 4 #' + 96G 4 MV)/ 

{[2GM - r) 3 r 5 ) (4.41) 

Higher powers of the d'Alembertian □ then lead to even more complicated expressions, with in- 
creasingly higher derivatives of a(r) and 9(r). Assuming a power law correction, as in Eqs. (4.28) 
and (4.29), with a = 0, as in Eq. (4.30), 

R = a MGa ( 2c a + c e {a - 1) ) r a ~ 3 + O{a 2 ) (4.42) 
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and then 

□ R -> a MGr a - 5 (2c a + c e (a - I)) a (a - 2) (a - 3) 

a 2 R -> a MGr a - 7 (2c a + c e (a- 1) ) a (a - 2) (a - 3) (a - 4) (a - 5) (4.43) 
and so on, and for general n — > 

a M G (2c a + c e (a-l)) "^J" " ~ ^ r"^' 3 (4.44) 

1 (2 — a) 

Therefore the only possible power solution for r 3> MG is a = 0, 2 ... i + 1, with c a and c# 
unconstrained to this order. 

4.3 Full effective field equations 

Next we examine the full effective field equations (as opposed to just their trace part) as in Eq. (4.15) 
with A = 0, 

area + AO) (JW-W) =° ("•«) 

valid for r / 0. If one denotes by 5G Ml/ = 5 [r^ v — \ Rj the lowest order variation (that is, of 
order cto) in the Einstein tensor over the ordinary vacuum solution G^ = 0, then one has 

again for r / 0. Here the covariant d'Alembertian operator 

□ = g^VvV* (4.47) 

acts on a second rank tensor, 

^ 7 uT a p = d u T a p — T^ u T X /3 - Tp u T a \ = I va p 

V M (V^T Q/3 ) = dplvap - ~ FapIuXP ~ ^^vaX (4-48) 

and would thus seem to require the calculation of as many as 1920 terms, of which many fortunately 
vanish by symmetry. In the static isotropic case the components of the Einstein tensor are given 
by 



A'B B B 

G -— — — 

B /2 r 2 A'B'r 2 B"r 2 A'r B'r 

^89 — A A T) < ?. T~H~F> + r> A TZ> ~ r> A 9 



4AB 2 AA 2 B 2AB 2A 2 2AB 

2 , 



G w = sin 2 9 G e0 (4.49) 
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Using again the expansion of Eqs. (4.25) and (4.26) one obtains for the tt, rr and 00 components 
of the Einstein tensor, to lowest order in a and 0, 

(2GM -r) a 1 



G tt ^ 



-0 + a + (-2GM + r) 0' 



— / r> /~t 11 /r \2 



Gee 



r(-2GM + r y 

{-GM+ r) (9-a+(2GM -r) {9' - a')) + r (-2 G M + rf 9" 



2{-2GM + ry 

G w = sin 2 9 Gee (4.50) 

After acting with □ on this expression one finds a rather complicated result. Here we will only list 
(□GV 

6BA' 3 2BA' 2 4B'A' 2 2BA' 6BA"A' B"A' 6B 
rA 5 + r 2 A 4 rA 4 ~ r 3 A 3 rA 4 + rA 3 + ^A 
6B 4B' 4B' BA" 2B'A" B" B" BA^ 

^4^2 r^yl V 3 A 2 V 2 A 3 tA 3 T^A ^2^2 j*A 3 ' 
To lowest order in cr(r) and 9{r) one finds the slightly simpler expressions 

(DG) tt = -(-2 G 2 M 2 9 + 2 G 2 M 2 a -4 G 3 M 3 9' + 2 G 2 M 2 r 6' - 28 G 3 M 3 a' + 38 G 2 M 2 r a' 
-16 GM r 2 a' + 2 r 3 a' + 24 G 3 M 3 r a" - 32 G 2 M 2 r 2 a" + 14 GM r 3 a" 
-2 r 4 a" - 8 G 3 M 3 r 2 + \2G 2 M 2 r 3 - 6 GM r 4 a® + r 5 c/ 3 ) ) / 
(r 6 (-2GM + r)) (4.52) 



(□G) rr = (6 G 2 M 2 9-4 GM r 9 - 6 G 2 M 2 a + 4 GM r a + 12 G 3 M 3 9' - 14 G 2 M 2 r 9' 

+4 GM r 2 0' - 12 G 3 M 3 a + 6 G 2 M 2 r a' + 4 GM r 2 a' - 2 r 3 a' + 8 G 3 M 3 r 0" 
-12 G 2 M 2 r 2 9" + 6 GM r 3 0" - r 4 0" + 4 G 2 M 2 r 2 a" - 4 GM r 3 a" + r 4 a" 
-8 G 3 M 3 r 2 (3) + 12 G 2 M 2 r 3 (3) - 6 GM r 4 (3) + r 5 (3) )/ 
(r 4 (-2GM + r) 3 ) (4.53) 



(□G) ee = (24 G 3 M 3 6»-36 G 2 M 2 r + 16 GM r 2 - 24 G 3 M 3 a + 36 G 2 M 2 r a - 16 GM r 2 a 

+48 G 4 M 4 0' - 96 G 3 M 3 r 0' + 68 G 2 M 2 r 2 9' - 16 GM r 3 9' + 16 G 4 M 4 cr' - 32 G 3 M 3 r a' 
+28 G 2 M 2 r 2 a - 16 GM r 3 a + 4 r 4 a' + 8 G 4 M 4 r 0" - 12 G 3 M 3 r 2 0" + 10 G 2 M 2 r 3 0" 
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-5 GM r 4 0" + r 5 6" - 24 G 4 M 4 r a" + 52 G 3 M 3 r 2 a" - 46 G 2 M 2 r 3 a" + 19 GM r 4 a" 
-3 r 5 a" - 24 G 4 M 4 r 2 # (3) + 44 G 3 M 3 r 3 # (3) - 30 G 2 M 2 r 4 # (3) + 9 GM r 5 # (3) - r 6 fl (3) 
+8 G 4 M 4 r 2 a (3) - 20 G 3 M 3 r 3 <r (3) + 18 G 2 M 2 r 4 <r (3) - 7 GM r 5 a (3) + r 6 a (3) 
+16 G 4 M 4 r 3 9^ - 32 G 3 M 3 r 4 #( 4 ) + 24 G 2 M 2 r 5 9^ - 8 GM r 6 #W + r 7 #( 4 ))/ 
(2r 3 (-2GM + r) 3 ) (4.54) 

with the if ip component proportional to the 66 component. If one again assumes that the corrections 
are given by a power, as in Eqs. (4.28) and (4.29), with a = (3, then one has to zeroth order 

G tt = a MGc a ar a - 3 

G rr = -a MG(c a + c e (a-l))r a - 3 

G ee = -\a Q MG{c c + c 6 {a-l)){a-l)r a - 1 (4.55) 

with the ipip component again proportional to the 66 component. Applying □ on the above Einstein 
tensor one then gets 

(□G) tt = a GMc a a(a-2){a-3)r a - 5 

(□G) rr = -a GM (c a + c e (a-l)) a(a-3)r a - 5 

{UG) ee = -^a GM (c a + cg(a-l)) a(a-3) 2 r a - 3 (4.56) 

again with the ipip component proportional to the 66 component. Applying □ again one obtains 

(□ 2 g)^ = a GMc a a(a-2)(a-3){a-4)(a-5)r a - 7 

(a 2 G) = -a GM (c, + C0 (a - 1)) a (a - 2) (a - 3) {a - 5) r a ~ 7 

\ / TV 

(u 2 G) = -\a GM (c a + c e {a-l)) a (a - 2) (a - 3) {a - 5) 2 r a ~ 5 (4.57) 
and so on, and for general n — > +^ 

(□ n G) tt -+ ^M, ^H fa -3 4 

(□™G) rr - - ao GM(c a + ce(a-l))- ^±t^±— r «-3-i 

(a - 1) (a - I) r(-a) 

i fa - 1 - l ) T(2 + i - a) 

(°"G), e - a GM(c CT + C ,(a-l))^ ^ » (4.58) 

Z (a-1) (^a--J r(-a) 

Inspection of the above results reveals a common factor l/r(— a), which would allow only integer 
powers a = 0, 1, 2 . . ., but the additional factor of l/(a — 1) excludes a = 1 from being a solution. 
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Even for a close to \jv (as expected on the basis of the non-relativistic expression of Eq. (3.24), as 
well as from Eq. (4.31)) v ~ 1/q — e only integer values a = 2, 3, 4 . . . are allowed. For the covariant 
divergences V^(n ra G)^ v one has 

V»(5G)^ r = (4.59) 

and, at the next order, 

V" (□ G% r = 2 a G 2 M 2 a (a - 3) (c a + c {\ - a)) r a ~ 7 (4.60) 

with the other components vanishing identically, and 

\7» (a 2 G) = 4 a G 2 M 2 a(a- 3) (a - 5) 2 (c CT + c<?(l - a)) r Q ~ 9 (4.61) 

again with the other components of the divergence vanishing identically. 

In general the problem of finding a complete general solution to the effective field equations by 
this method lies in the difficulty of computing arbitrarily high powers of □ on general functions 
such as a(r) and 6(r), which eventually involve a large number of derivatives. Assuming for these 
functions a power law dependence on r simplifies the problem considerably, but also restricts the 
kind of solutions that one is likely to find. More specifically, if the solution involves (say for small 
r, but still with r 3> 2MG) a term of the type r a lnmr, as in Eqs. (3.21), (5.35) and (5.38) for 
v — ► 1/3, then this method will have to be dealt with very carefully. This is presumably the reason 
why in some of the T-function coefficients encountered here one finds a power solution (in fact 
a = 3) for v close to a third, but one gets indeterminate expression if one sets exactly a = \ jv = 3. 



5 The Quantum Vacuum as a Fluid 



The discussion of Sec. 3 suggests that the quantum correction due to the running of G can 
be described, at least in the non-relativistic limit of Eq. (2.2) as applied to Poisson's equation, 
in terms of a vacuum energy density p m (r), distributed around the static source of strength M 
in accordance with the result of Eqs. (3.17) and Eq. (3.19). These expressions, in turn, can be 
obtained by Fourier transforming back to real space the original result for G(k 2 ) of Eq. (2.2). 

Furthermore, it was shown in Sec. 4 (and was discovered in [12] as well, see for example Eq. (7.8) 
later in this paper) that a manifestly covariant implementation of the running of G, via the G(D) 
given in Eqs. (4.2) and (4.7), will induce a non- vanishing effective pressure term. This result can 
be seen clearly, in the case of the static isotropic metric, for example from the result of Eq. (4.35). 
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We will therefore, in this section, consider a relativistic perfect fluid, with energy-momentum 
tensor 

T/m, = [p + p] u^u v + g^p (5.1) 
which in the static isotropic case reduces to Eq. (4.13), 

= diag [ B (r) p(r), A{r)p(r), r 2 p(r), r 2 sin 2 9 p(r)] (5.2) 

and gives a trace T = 3p — p. 

The tt, rr and 99 components of the field equations then read 

-ABM + - + ^ = SnGB(r)p(r) (5.3) 

v ' rA(r) 2 r 2 A(r) r 2 u w y J 

XA(r) - ^ + ^ + 1 = 8,GiM P (r) (5.4) 

g(g I Ar 2 A'(r)i?'(r)r 2 £"(r)r 2 A'(r)r ^(r)r o r 2 frl f( . ^ 

4A(r)5(r) 2 4,4(r) 2 5(r) 2A(r)5(r) 2.4(r) 2 + 2i(r)B(r) PW 1 j 

with the 9999 component equal to sin 2 9 times the 99 component. 
Energy conservation V Ij^ = implies 

[pW + pWI^+p'W = (5.6) 

and forces a definite relationship between B(r), p{r) and p(r). The three field equations and the 
energy conservation equation are, as usual, not independent, because of the Bianchi identity. 

It seems reasonable to attempt to solve the above equations (usually considered in the context 
of relativistic stellar structure [67]) with the density p(r) given by the p m (r) of Eqs. (3.17), (3.18) 
and (3.19). 

This of course raises the question of how the relativistic pressure p(r) should be chosen, an issue 
that the non-relativistic calculation did not have to address. We will argue below that covariant 
energy conservation completely determines the pressure in the static case, leading to consistent 
equations and solutions (note that in particular it would not be consistent to take p(r) = 0). 

Since the function B(r) drops out of the tt field equation, the latter can be integrated immedi- 
ately, giving 

A{r)- 1 = l + £l_^ r 2_^ f d xx 2 p(x) (5.7) 
r 3 r Jo 
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which suggests the introduction of a function m(r) 

rr 

m(r) = 47r / dxx 2 p(x) (5-8) 

JO 

It also seems natural in our case to identify C\ = —2MG, which of course corresponds to the 
solution with ao = (p = p = 0) (equivalently, the point source at the origin of strength M could 
be included as an additional 5-function contribution to p(r)). 
Next, the rr field equation can be solved for B(r), 



B(r) = exp L - f dy ^^^11} (5 . 9) 



l + A(y)(Xy 2 -8 7TGy 2 p(y) - 1) 

/r o y 

with the constant C2 again determined by the requirement that the above expression for B(r) reduce 
to the standard Schwarzschild solution for ao = (p = p = 0), giving C2 = ln(l — 2MG/ro — Atq/3). 
The last task left therefore is the determination of the pressure p(r). 

Using the rr field equation, B'(r)/B(r) can be expressed in term of A(r) in the energy conser- 
vation equation, which results in 

2rp'(r)- [l + A(r)(Ar 2 -87rGr 2 p(r)-l)j (p(r)+p(r)) =0 (5.10) 

Inserting the explicit expression for A(r), from Eq. (5.7), one obtains 

(8irG r 3 p(r) + 2MG - |Ar 3 + 8vrG [l dxx 2 p(x)) (p(r) + pir)) 

p'(r) + A 1 °- J = (5.11) 

2 r (r - 2MG - | r 3 - 8itG J r dx x 2 p(x)J 

which is usually referred to as the equation of hydrostatic equilibrium. Prom now on we will focus 

only the case A = 0. Then 

(SttG r 3 p(r) + 2MG + 8nG f£ dxx 2 p(x)) (p(r) + p{r)) 
P[r)+ 2r(r - 2MG - 8irG Jq dx x 2 p(x)) [ ' 

The last equation, a non-linear differential equation for p(r), can be solved to give the desired 
solution p(r), which then, by equation Eq. (5.9), determines the remaining function B(r). 

In our case though it will be sufficient to solve the above equation for small ao, where ao (see 
Eq. (2.2) and Eq. (3.17)) is the dimensionless parameter which, when set to zero, makes the solution 
revert back to the classical one. 

It will also be convenient to pull out of A(r) and B(r) the Schwarzschild solution part, by 
introducing the small corrections a(r) and 0(r) (already defined before in Eqs. (4.25) and (4.26)), 
both of which are expected to be proportional to the parameter ao- One has 

a(r) = -8ttG rdxx 2 p(x) = -2m(r)G (5.13) 
Jo 
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and 

* M = » P {" + £ rf %-2MG + -8^1^(x) } + 2MG - r (5 ' 14) 
Again, the integration constant C2 needs to be chosen here so that the normal Schwarzschild solution 

is recovered for p = p = 0. 

To order ao the resulting equation for p(r), from Eq. (5.12), is 

Note that in regions where p(r) is slowly varying, p'(r) ~ 0, one has p ~ — p, i.e. the fluid 
contribution is acting like a cosmological constant term with a(r) ~ 6{r) ~ — (p/3)r 3 . 
The last differential equation can then be solved for p(r), 

where the constant of integration has to be chosen so that when p(r) = (no quantum correction) 
one has p(r) = as well. Because of the singularity in the integrand at r = 2MG, we will take the 
lower limit in the integral to be rg = 2MG + e, with e — > 0. 

To proceed further, one needs the explicit form for p m (r), which was given in Eqs. (3.17), (3.18) 
and (3.19). The required integrands involve for general v the modified Bessel function K n (x), and 
can be therefore a bit complicated. But in some special cases the general form of the density p m 
of Eq. (3.17) 

Pm(r) = ^-c u a Mm 3 (mr)'^'^ ifi (3 .u(mr) (5.17) 
reduces to a relatively simple expression, which we will list here. For v = 1 one has 

1 r, 1 

Pm(r) = —^a Mm 6 Ki(mr) (5.18) 

2tt z mr 



whereas for v = 1/2 one has 

and for v = 1/3 

and finally for v = 1/4 



p m (r) = — a Mm 3 — e~ mr (5.19) 
47r mr 



p m (r) = ^a Mm 3 K (mr) (5.20) 



p m (r) = J-a Mm 3 e~ mr (5.21) 

87T 



Note that p m (r) diverges at small r for v > 1/3 

Here we will limit our investigation to the small r (mr <^ 1) and large r (mr 1) behavior. 
Since m = l/£ is very small, the first limit appears to be of greater physical interest. 
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5.1 Small r limit 

For small r the density p m (r) has the following behavior (see Eq. (3.17)), 

Pm(r) ~A r^- 3 (5.22) 

r— > 

for v > 1/3, with 

ci~ c i sec ( ~ ) i 

A) = jLJi a Mmv = — , 2vJ \ a Mm- (5.23) 

87r 47rr(i-lJ 

where the dimensionless positive constant c& is determined from the small x behavior of the modified 
Bessel function K n (x), 

wi ,\ , . 22 7r sec (t?-) i o i o , 

x^i^-VKif. M (x) ~ 7 r^x-" 3 = c fc x-~ 3 5.24 

1 \2v 2 J 

valid for v > 1/3, and is given in Eq. (3.18). For u < 1/3 p m (r) ~ const, clq Mm 3 , independent 
of r. For v = 1/3 the expression for p m {r) in Eq. (5.20) should be used instead. 
Therefore in this limit, with | < v < 1, one has 

m(r) ~ AttvAqt~ (5.25) 

and, from the definition of <r(r), 

cr(r) ~ -2m(r)G = -%-kvGAqt^ (5.26) 

and finally 

2 ikf (j i i 

A~ l (r) = 1 2a MGc s m«rH + ... (5.27) 



with the constant c s = u\ sec (^) |/r(£ — 1). For f = 1/3 the last contribution is indistinguishable 



from a cosmological constant term — ^r 2 , except for the fact that the coefficient here is quite 



different, being proportional to ~ aoMGm 3 . 

To determine the pressure, we suppose that it as well has a power dependence on r in the regime 
under consideration, p m (r) = c p A r 7 , where c p is a numerical constant, and then substitute p m (r) 
into the pressure equation Eq. (5.15). This gives, past the horizon r S> 2MG, 

(2 7 - l)c p MGr"<- 1 - c pl r^ - MGr l ' u ~ 4 ~0 (5.28) 

giving the same power 7 = 1/v — 3 as for p(r), c p = — 1 and surprisingly also 7 = 0, implying that in 
this regime only v = 1/3 gives a consistent solution. Again, the resulting correction is quite similar 
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to what one would expect from a cosmological term, with an effective A m /3 ~ ^ v a^M Gmv . 
One then has for v near 1/3 

Pm (r) = A c p r^~ 3 + ... (5.29) 

and thus from Eq. (5.15) 

B(r) = 1 - ^£J1 _ 2a MGc s mv r*- 1 + ... (5.30) 
r 

Both the result for A(r) in Eq. (5.27), and the above result for B{r) are, for r » 2M G, consistent 
with a gradual slow increase in G in accordance with the formula 



G 



G(r) = G (l + a c s mv r» + . . .) (5.31) 



We note here that both expressions for A(r) and B(r) have some similarities with the approximate 
non-relativistic (Poisson equation) result of Eq. (3.24), with the correction proportional to ao 
agreeing roughly in magnitude (but not in sign). 

The case v = 1/3 requires a special treatment, since the coefficient c& in Eq. (5.24) diverges as 
v — > 1/3. Starting from the expression for p m (r) for v = 1/3 in Eq. (5.20), 

1 



p m (r) = aoMm 6 K (mr) (5.32) 



one has for small r 



Pm{r) = -^Mm 6 I In— + 7 )+ ... (5.33) 



and therefore from Eq. (5.14), 



3vr 

and consequently 



, , 4a MGm 3 , . . . . 
a(r) = r m(mr) + ... (5.34) 



,_ i- , , 2MG AclqM Gm 6 2l . , . or , 

A 1 (r) = 1 1 r 2 In (mr) + ... (5.35) 

r 37r 

From Eq. (5.15) one can then obtain an expression for the pressure p m {r), and one finds 

w 3l / \ a Mm 3 log ( r + n/T^G-MG 3 ,, 3 

. , a Q Mm A log(mr) u 6 I 1 V r / a Mm d a Mm 6 c 3 

P™{r) = r-o ; ; = ~+ o + 



2;r 2 2-^ 1 ^ ^ 2-^ 1 ■ n,<: 



(5.36) 

where C3 is again an integration constant. Here we will be content with the r 3> 2MG limit of the 
above expression, which we shall write therefore as 

Pm(r) = pMm 3 hH + ... (5.37) 
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After performing the required r integral in Eq. (5.14), and evaluating the resulting expression in 
the limit r 3> 2MG, one obtains an expression for 9(r), and consequently from it 

2MG 4a MGm 3 2n , 

B(r) = 1 + — — r 2 In (mr) + ... (5.38) 

r 3 7r 

The expressions for A(r) and i?(r) are, for r > 2 M G, consistent with a gradual slow increase in 
G in accordance with the formula 

G - G(r) = G (l + ^m 3 r 3 In + ..) (5.39) 

and therefore consistent as well with the original result of Eqs. (2.1) or (2.2), namely that the 
classical laboratory value of G is obtained for r <C £. In fact it is reassuring that the renormalization 
properties of G(r) as inferred from A(r) are the same as what one finds from B(r). Note that the 
correct relativistic small r correction of Eq. (5.39) agrees roughly in magnitude (but not in sign) 
with the approximate non-relativistic, Poisson equation result of Eq. (3.21). 

One further notices some similarities, as well as some rather substantial differences, with the 
corresponding QED result of Eq. (2.67). In the gravity case, the correction vanishes as r goes to zero: 
in this limit one is probing the bare mass, unencumbered by its virtual graviton vacuum polarization 
cloud. On the other hand, in the QED case, as one approaches the source one is probing the bare 
charge, unscreened by the electron's vacuum polarization cloud, and whose magnitude diverges 
logarithmically for small r. 

It should be recalled here that neither function A(r) or B(r) are directly related to the relativistic 
potential for particle orbits, which is given instead by the combination 



Veff(r) 



2A(r) 



• 2 1 +i 



.,2 



B(r) 



(5.40) 



where I is proportional to the orbital angular momentum of the test particle [68] . 

Furthermore, from the metric of Eqs. (5.27) and (5.30) one finds for v —> 1/3 the following 
results for the curvature invariants 

R 2 = 1024^G 2 7r 2 



R^Ri™ = 256 AiG 



R^R^ = 16G 2 (^^ + ^) (5.41) 

which are non-singular at r = 2 M G, and again consistent with an effective mass density around 
the source m(r) oc r 3 . 
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5.2 Large r limit 

For large r one has instead, from Eq. (3.17) for p m (r), 

Pm {r) ~ A r^~ 2 e- mr (5.42) 

r— > oo 

with = 1/V 1287T cv ao M m 1+ ^ . In the same limit, the integration constants is chosen so that 
the solution for A(r) and B(r) at large r corresponds to a mass M 1 = (l + an)M (see the expression 
for the integrated density in Eq. (3.19)), or equivalently 

a( r ) ~ 9(r) ~ -2a MG (5.43) 

r — > oo 

On then recovers a result similar to the non-relativistic expression of Eqs. (3.7), (3.8) and (3.12), 
with G(r) approaching the constant value = (1 + ao)G, up to exponentially small corrections 
in mr at large r. 

In conclusion, it appears that a solution to relativistic static isotropic problem of the running 
gravitational constant can be found, provided that the exponent v in either Eq. (2.2) or Eq. (4.8) 
is close to one third. This last result seems to be linked with the fact that the running coupling 
term acts in some way like a local cosmological constant term, for which the r dependence of the 
vacuum solution for small r is fixed by the nature of the Schwarzschild solution with a cosmological 
constant term 13 . 



6 Distortion of the Gravitational Wave Spectrum 



A scale-dependent gravitational constant G(k 2 ) will cause slight distortions in the spectrum of 
gravitational radiation at extremely low frequencies, to some extent irrespective of the nature of 
the perturbations that cause them. From the field equations with A = 

iV -\dp,vR = 8ttGT^ (6.1) 
one obtains in the weak field limit with harmonic gauge condition 

□ V = SttGT^ (6.2) 

13 In d > 4 dimensions the Schwarzschild solution to Einstein gravity with a cosmological term is [69] A _1 (r) = 
B(r) = 1 - 2MGc d r 3 ~ d - (d _ 2 ^ d _ 1) r 2 , with c d = 47iT(^-)/(d - 2)tt^~ , which would suggest, in analogy with the 
results for d = 4 given in this section, that in d > 4 dimensions only v = l/(d — 1) is possible. This last result would 
also be in agreement with the exact value v = found at d — oo [15] 
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with as usual 

%v = - \r,^ u T x x (6.3) 

Density perturbations <5p(x, t) will enter the r.h.s. of the field equations and give rise to gravitational 
waves with Fourier components 

V(*0 = -87rGpT^(p,p)(fc) (6.4) 
giving for the power spectrum of transverse traceless (gravitational wave) modes 

P TT (k 2 ) ~ k 3 \h TT (k)\ 2 = (8Tr) 2 G 2 ±\T(p,p)(k)\ 2 (6.5) 
A scale dependent gravitational constant, with variation in accordance with Eq. (2.2), 

G -► G(k 2 ) (6.6) 
would affect the spectrum of very long wavelength modes via 



P TT (k 2 ) ~ fc 3 |V(*)| 2 = (8^) 2 G 2 (fc 2 )i|f^(p,p)| 2 (6.7) 



Specifically, according to the expression in Eq. (2.2) for the running of the gravitational constant, 

G(k 2 ) ( m 2 \& 

+ - (6 - 8) 

one has for the tensor power spectrum 



P TT {k 2 ) ~ k 3 \h^(k)\ 2 = (8n) 2 G 2 ^- 



/ 771 

1 + ao [WT 



i 

2 \ 2^ 



2 

2 



T^fop)! 2 (6.9) 



with the expression in square brackets varying perhaps by as much as an order of magnitude from 
short wavelengths k S> l/£, to very long wavelengths A; ~ l/£. 

7 Quantum Cosmology - An Addendum 

In this section we will discuss briefly what modifications are expected when one uses Eq. (2.2) 
instead of Eq. (2.1) in the effective field equations. In [12] cosmological solutions within the 
Friedmann- Robertson- Walker (FRW) framework were discussed, starting from the quantum ef- 
fective field equations of Eq. (2.52), 

iV ~ h^R + Xg^ = 8vrG (l + A(a)) (7.1) 
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with A(0) defined in either Eq. (4.6) or Eq. (4.7), and applied to the standard Robertson- Walker 
metric 



-dt 2 + a 2 {t) 



dr z 



+ r z 



(d6 2 + sm 2 6d^ 



(7.2) 



, 1 — k r 2 

It should be noted that there are two quantum contributions to this set of equations. The first one 
arises because of the presence of a non- vanishing cosmological constant A ~ l/£ 2 , as in Eq. (2.48), 
originating in the non-perturbative vacuum condensate of the curvature. As in the case of standard 
FRW cosmology, this is the dominant contributions at large times t, and gives an exponential 
expansion of the scale factor. 

The second contribution arises because of the running of G for t <C £ in the effective field 
equations, 



+ ... 



(7.3) 



G(a) = G (1 + A(a)) = G 

with v ~ 1/3 and ao a calculable coefficient of order one (see Eqs. (2.1) and (2.2)). 

In the simplest case, namely for a universe filled with non-relativistic matter (p=0), the effective 
Friedmann equations then have the following appearance [12] 



k a 2 {t) 
a 2 (t) + a 2 (t) 

for the tt field equation, and 
k 



8vrG(t) 
3 

8vrG 



1 



M + 3e 

1 + c C (t/0 1/u + 



p{t) + W 2 



a 2 ^ 2a(t) 



8vrG 



c 5 (t/0 



i/i/ 



+ 



p(t) + 



(7.4) 



(7.5) 



a 2 (t) a 2 (t) ' a(t) 

for the rr field equation. The running of G appropriate for the RW metric, and appearing explicitly 
in the first equation, is described by 



G(t) = G 



1 + « I e 



+ ... 



(7.6) 



(with eg or the same order as ao of Eq. (2.1) [12]). Note that the running of G(t) induces as well 
an effective pressure term in the second (rr) equation. 14 
One can therefore talk about an effective density 

G(t) 



Peff(t) 



G 



Pit) 



(7.7) 



14 We wish to emphasize that we are not talking here about models with a time-dependent value of G. Thus, for 
example, the value of G ~ G c at laboratory scales should be taken to be constant throughout most of the evolution 
of the universe. 
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and an effective pressure 



Peff(t) 



1 /G(t) 



- l) Pit) 



(7.8) 



3 V G 



with p e ff(t)/p e ff(t) = \{G{t) — G) / G{t) 15 . Within the FRW framework, the gravitational vacuum 
polarization term behaves therefore in some ways (but not all) like a positive pressure term, with 
pit) = u> p(t) and uj = 1/3, which is therefore characteristic of radiation. One could therefore 
visualize the gravitational vacuum polarization contribution as behaving like ordinary radiation, 
in the form of a dilute virtual graviton gas: a radiative fluid with an equation of state p = It 
should be emphasized though that the relationship between density p(t) and scale factor a(t) is 
very different from the classical case. 

The running of G(t) in the above equations follows directly from the basic result of Eq. (2.1) 
(with the dimensionless constant proportional to do, with a numerical coefficient of order one 
given in magnitude in [12]), but transcribed, by explicitly computing the action of the covariant 
d'Alembertian □ = g^ v V M Vj, on T^, for the RW metric. In other words, following the more or 
less unambiguously defined sequence G(k 2 ) — > (?(□) — > G(t). At the same time, the discussion of 
Sec. 1 underscores the fact that for large times t>^ the form of Eq. (2.1), and therefore Eq. (7.6), 
is no longer appropriate, due to the spurious infrared divergence of Eq. (2.1) at small k 2 . Indeed 
from Eq. (2.2), the infrared regulated version of the above expression should read instead 



with £ = m^ 1 the (tiny) infrared cutoff. Of course it reduces to the expression in Eq. (7.6) in 
the limit of small times t, but for very large times f > ( the gravitational coupling, instead of 
unphysically diverging, approaches a constant, finite value Goo = (1 + ao + . . .) G c , independent of 
£. The modification of Eq. (7.9) should apply whenever one considers times for which t <C £ is not 



of academic interest, and was therefore not discussed much in [12]. 

It should be noted that the effective Friedman equations of Eqs. (7.4) and (7.5) also bear 
a superficial degree of resemblance to what might be obtained in some scalar-tensor theories of 
gravity, where the gravitational Lagrangian is postulated to be some singular function of the scalar 

curvature [70, 71]. Indeed in the FRW case one has, for the scalar curvature in terms of the scale 

Strictly speaking, the above results can only be proven if one assumes that the pressure's time dependence is 
given by a power law, as discussed in detail in [12]. In the more general case, the solution of the above equations for 
various choices of £ and ao has to be done numerically. 




(7.9) 



valid. But since £ ~ is of the order the size of the visible universe, the latter regime is largely 
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factor, 

R = 6(k + a 2 (t) + a(*)a(t)) /a 2 (t) (7.10) 
and for /c = and a(i) ~ t a one has 

ij = 6a(2 ;~ 1} (7.ii) 

which suggests that the quantum correction in Eq. (7.4) is, at this level, nearly indistinguishable 
from an inverse curvature term of the type (£ 2 .R) -1 / 2 ^, or 1/(1 + £ 2 R) 1 / 2v if one uses the infrared 
regulated version. The former would then correspond the to an effective gravitational action 

^' a wsh^( R + -^- tx ) (7 ' 12) 

with / a numerical constant of order one, and A ~ l/£ 2 . But this superficial resemblance is seen 
here more as an artifact, due to the particularly simple form of the RW metric, with the coincidence 
of several curvature invariants not expected to be true in general. 



8 Conclusions 



In this paper we have examined a number of basic issues connected with the renormalization 
group running of the gravitational coupling. The scope of this paper was to explore the overall 
consistency of the picture obtained from the lattice, by considering a number of basic issues, one 
of which is the analogy, or contrast, with a much better understood class of theories, namely QED 
on the one hand, and non-abelian gauge theories and QCD on the other. 

The starting point for our discussion of the renormalization group running of G (Sec. 2) is 
Eq. (2.1) (valid at short distances k 3> m, or, equivalently r <C £), and its improved infrared 
regulated version of Eq. (2.2). The scale dependence for G obtained from the lattice is remarkably 
similar to the result of the 2 + e expansion in the continuum, as in Eq. (2.59), with two important 
differences: only the strong coupling phase G > G c is physical, and for the exponent one has 
v ~ 1/3 in four dimensions. The similarity between the two results in part also originates from 
the fact that in both cases the renormalization group properties of G are inferred (implicitly, in 
the 2 + e case) from the requirement that the non-perturbative scale of Eq. (2.10) be treated as an 
invariant. 

Inspection of the quantum gravitational functional integral Z of Eq. (2.20) reveals that its 
singular part can only depend of the dimensionless combination AqG 2 , up to an overall factor 
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which cannot affect the non-trivial scaling behavior around the fixed point, since it is analytic in 
the couplings. This then leaves the question open of which coupling(s) run and which ones do not. 

The answer in our opinion is possibly quite simple, and is perhaps best inferred from the nature 
of the Wilson loop of Eq. (2.47): the appropriate renormalization scheme for quantum gravity is 
one in which G runs with scale according to the prediction Eq. (2.2), and the scaled cosmological 
constant A is kept fixed, as in Eqs. (2.48) and (2.51). Since the scale £ is related to the observable 
curvature at large scales, it is an almost inescapable conclusion of these arguments that it must 
be macroscopic. Furthermore, it is genuinely non-perturbative and non-analytic in G, as seen for 
example from Eq. (2.54), and represents the effects of the gravitational vacuum condensate which 
makes its appearance in the strongly coupled phase G > G c . 

Another aspect we have investigated in this paper is the nature of the quantum corrections to the 
gravitational potential <p(r) in real space, arising from the scale dependence of Newton's constant 
G. The running is originally formulated in momentum space (see Eq. (2.2)), since it originates in 
the momentum dependence of G as it arises on the lattice, or in the equivalent renormalization 
group equations, Eqs. (2.9) or (2.17). The solution <j>(r) to the non-relativistic Poisson equation for 
a point source is given in Eq. (3.21) of Sec. 3 for various values of the exponent v. The solution 
is obtained by first computing the effective vacuum polarization density p m (r) of Eq. (3.17), and 
then using it as a source term in Poisson's equation. Already in the non-relativistic case, the value 
v = 1/3 appears to stand out, since it leads to logarithmic corrections at short distances r <C £. 

A relativistic generalization of the previous results was worked out in Sees. 4 and 5. First it 
was shown that the scale dependence of G can be consistently embedded in a relativistic covariant 
framework using the d'Alembertian □ operator, leading to a set of nonlocal effective field equations, 
Eq. (4.8). The consequences can then be worked out in some detail for the static isotropic metric 
(Sec. 4), at least in a regime where 2MG <C r <C £, and under the assumption of a power law 
correction (otherwise the problem becomes close to intractable). One then finds that the structure 
of the leading quantum correction severely restricts the possible values for the exponent v, in the 
sense that no consistent solution to the effective non-local field equations, incorporating the running 
of G, can be found unless v~ x is an integer. 

A somewhat different approach to the solution of the static isotropic metric was then discussed 
in Sec. 5, in terms of the effective vacuum density of Eq. (3.17), and a vacuum pressure chosen 
so as to satisfy a covariant energy conservation for the vacuum polarization contribution. The 
main result is the derivation from the relativistic field equations of an expression for the metric 
coefficients A(r) and B(r), given in Eqs. (5.35) and (5.38). For v = 1/3 it implies for the running 
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of G in the region IMG < r <C £ the result of Eq. (5.39), 

G(r) = G (l + ^m 3 r 3 In -i^ + . . (8.1) 
V 3 7r m^r / 

indicating therefore a gradual, very slow increase in G from the "laboratory" value G = G(r = 0). 
For the actual values of the parameters appearing in the above expression one expects that m is 
related to the curvature on the largest scales, m~ 1 = £ ~ 10 28 cm, and that a ~ O(10). From 
the nature of the solution for A(r) and B(r) one finds again that unless the exponent v is close to 
1/3, a consistent solution of the field equations cannot be found. Note that for very large r 3> £ 
the growth in G(r) saturates and the value G^ = (1 + ao)G is obtained, in accordance with the 
original formula of Eq. (2.2) for k 2 ~ 0. A natural comparison is with the QED result of Eq. (2.67). 

At the end of the paper we have added some remarks on the solution of the gravitational wave 
equation with a running G. We find that a running Newton's constant will slightly distort the gra- 
vitational wave spectrum at very long wavelengths (Sec. 6), according to Eq. (6.9). Regarding the 
problem of finding solutions of the effective non-local field equations in a cosmological context [12], 
wherein quantum corrections to the Robertson- Walker metric and the basic Friedman equations 
(Eqs. (7.4) and (7.5)) are worked out, we have discussed some of the simplest and more plausible 
scenarios for the growth (or lack thereof) of the coupling at very large distances, past the de Sitter 
horizon. 
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